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Abstract: Using the AdS/CFT correspondence, we study the anisotropic charge trans- 
port properties of both supersymmetric and non-supersymmetric matter fields on (2+1)- 
i dimensional defects coupled to a (3+l)-dimensional J\f = 4 SYM "heat bath". 

qj We focus on the case of a finite external background magnetic field, finite net charge density 

and finite mass. At high frequencies, we discover a spectrum of quasiparticle resonances due 
to the magnetic field and finite density and at small frequencies, we perform a Drude-like 
expansion around the DC limit. Both of these regimes display many generic features and 
some features that we attribute to strong coupling, such as a minimum DC conductivity 
and an unusual behavior of the "cyclotron" and plasmon frequencies, which become related 
to the resonances found in the conformal case in an earlier paper. We further study the 
hydrodynamic regime and the relaxation properties, from which the system displays a set 
of different possible transitions to the collisionless regime. The mass dependence can be 
cast in two regimes: a generic relativistic behavior dominated by the UV and a non-linear 
J> hydrodynamic behavior dominated by the IR. In the massless case, we furthermore extend 

earlier results from the literature to find an interesting selfduality under a transformation 
^ of the conductivity and the exchange of density and magnetic field. 
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1. Introduction 



The AdS/CFT correspondence [1, 2, 3] is a very important tool to study strongly coupled 
field teories using both string theory based and other gravitational setups. The first, and 
most common, example is the duality between a (black) stack of N c D3 branes generating 
an AdS§ x S* 5 geometry in the decoupling limit and a thermal J\f = 4 SU(N C ) super- 
Yang-Mills theory on its boundary Of particular importance to experiment has been the 
introduction of fields transforming in the fundamental representation of the SU (N c ) by 
introducing a small number of Nf "probe" D7-branes into the background, covering all 
of the AdS directions [4, 5]. The U(Nf ) symmetry of the stack of probe branes then 
corresponds to the flavor symmetry. Obviously this setup is still significantly different 
from the standard model QCD, but studying the physical properties of both this model 
and the T-dual D4-D6 setup have received great interest [6, 7]. Particular activity has been 
related to the thermodynamics and the phase structure [8] and to the "meson spectrum" 
[9]. The hope for experimental matchings is that some results obtained from AdS/CFT 
may be sufficiently generic, such that they also apply to QCD. 

Recently there has also been great effort on applying the AdS/CFT correspondence 
to condensed matter physics. One particular motivation for why this may be successful 
is the fact that while there is only one QCD, there are on the one hand many different 
strongly coupled effective field theories in condensed matter physics and on the other hand 
there exist a large number AdS string vacua. Hence thre is a significant potential for 
close matchings of gravitational duals to experimentally relevant CFTs. Starting with 
the study of transport properties of 2+1 dimensional field theories dual to an AdS^ x S" 6 
geometry obtained from an M2-brane setup in [10], there has been a huge amount of 
activity studying 2+1 dimensional systems. Those systems display a variety of interesting 
properties, such as Hall conductivity [11, 12], superconductivity and superfluidity [13]. 
While most constructions are not based on string theory, thre are some string theory based 
constructions that ensure that the systems are pathology-free. A good example for how 
realistic those constructions are is e.g. the reproduction of the quantization of the magnetic 
flux in superconductors [14] or the Nernst effect [15]. Another interesting example of recent 
progress is the contruchtion of duals of non-relativistic CFTs [16]. A yet not satisfactorily 
addressed question is what the implications of the Fermi-Dirac distribution is in purely 
fermionic systems [17]. 

The most common example of how systems with conformal symmetry arise in con- 
densed matter physics is the quantum critical phase. This phase arises in the context of 
a phase transition at zero temperature, a so-called "quantum critical phase transition" at 
a "quantum critical point". At this point, the system displays sclae-invariant behavior as 
also in other phase transitions, but in contrast to phase transitions at finite temperature, 
it extends into a whole region in the phase diagram that may be described by a confor- 
mal field theory, the so-called quantum critical phase; shown in fig. 1. As an example, 
the conductivity in the quantum critical phase is thought to be controlled by a univer- 
sal function S that depends on the ratio oj/T and some dimensional temperature scaling, 
<j(lu) = <5 2 (T/c) rf ~ 2 X(o;/T) for some microscopically determined velocity c. Obtaining £ 
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however is a difficult task and can only be done in certain limits. For instance for large 
frequencies, one expects a ~ u d ~ 2 [18, 19]. 

All of the above-mentioned systems have in common that they are described by 2+1 
dimensional field theories. In our 3+1 dimensional world however, all 2+1 dimensional 
systems are strictly speaking defects. In some cases this fact may be less relevant and in 
other cases more relevant. Hence, it is interesting to study the physics of a 2+1 dimensional 
defect in order to explore what difference there is to purely 2+1 dimensional systems. 
Defect field theories are basically field theories in which matter that is confined to some 
hypersurface interacts via a field theory in the bulk. While there is some review literature 
in a soft condensed matter context, related to aspects like the statistical mechanics of 
crystal defects in the context of melting behaviors, defects in polymers or flux tubes in 
superconductors (for a review see [20]) there seems to be not much review literature related 
to the defects and their aspects that we are interested in. Hence, a motivated guess may 
be that they have many properties in common with surfaces, which have been studied 
extensively. To illustrate their properties, we can look in fig. 1 at a generic surface phase 
diagram of some system described by a bulk coupling Jj, and a surface coupling J s - for 
a review on the subject see [21]. There we see that over most of the parameter range the 
surface and the bulk are in the same phase and display a simultaneous "ordinary" phase 
transition. As we tune the surface coupling beyond a "special point" , which is some critical 
multiple of the bulk coupling, the phase transitions on the surface and in the bulk separate 
into a surface phase transition and an "extraordinary" phase transition in which there is 
a phase transition only in the bulk. It is obvious from the ratio J s /Jb in this regime that 
the ordered phase on the surface extends to higher temperatures than the ordered phase 
in the bulk. However it is quite interesting that this splitting of phase transitions typically 
occurs as J s becomes greater than J and hence there is no "mirror symmetric" version of 
this plot. 
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Figure 1: Left: A generic quantum critical phase transition and quantum critical phase. Right: 
A generic surface phase diagram. 



In this paper, we will study the transport properties of the matter along the defect 
in order to infer on its physical properties. In that sense we will be put sometimes in the 
perspective of an experimentalist, trying to interpret our results. In order to identify the 
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characteristics that are due to having a defect rather than just a plain 2+1 dimensional 
field theory, a major theme in this research is the inclusion of a parameter that is related 
to a difference in the level of the gauge group N c — ► 5N C between both sides of the defect. 
The particular focus is then to extend earlier results obtained in the purely conformal case 
to the case of finite background parameters. 

Also, we will try to link the properties obtained in different regimes rather than study- 
ing only one particular limit or one particular effect. In this spirit, we will study all 
frequency and wavenumber scales, from the DC limit and the hydrodynamic regime (the 
"collision-dominated" regime at small frequencies and small wavenumbers) up to the quasi- 
particle limit (the "collisionless" regime at large frequencies) - and we will demonstrate in 
several cases how a length scale obtained in one regime, direction or context will govern 
some properties also in another regime. Furthermore, we will be interested how effects 
that we are familiar with from ordinary weakly coupled free electron gas type physics will 
manifest themselves in this strongly coupled system, as we turn on various "condensed 
matter" parameters, i.e. the net baryon number density, background magnetic field and 
"quark" mass. Certainly, since we are working at finite temperatures, there will aways be a 
finite total "quark" density, and the net baryon number density in some sense corresponds 
to the difference of the number of "quarks" and "anti-quarks" . While those quantities may 
move us away from the quantum critical point in the defect field theory, the bulk theory 
will still remain J\f = 4 SYM. In terms of the phase diagram of surface phase transitions, 
this would move us along the direction of surface coupling towards the extraordinary phase 
transition. 

One interest is furthermore to explore what happens in this defect setup to the result 
of the constant conductivity due to the electromagnetic duality that was studied in [10], 
in particular at finite background quantities. 

Our defect CFT is realized by inserting Nf probe D5- or D7-branes into the background 
of a black D3-brane. In either construction, the difference 5N C in the level of the gauge 
group of the 3+1 SYM will be introduced by an additional flux on the probe brane in 
the compact sphere. The defect CFT constructed with the D5-branes is certainly well 
known [4, 22, 23]. Certain aspects of the D7-brane construction have also been studied 
previously [24, 25] but we should note that the internal flux introduced here is essential to 
remove an instability that would otherwise appear in this construction. The finite magnetic 
field and net density are introduced using the well-known duals of a magnetic field and 
an electric field, respectively, in the world volume of the probe brane. The finite quark 
mass will be obtained by a deformation of the embedding in the compact sphere in the 
same fashion in which it was done in the duals for 3+1 dimensional QCD-like systems 
[8, 26, 27, 7]. To obtain the transport properties, we will use linear response theory to 
study the conductivity on the defect at finite frequency and temperature and at finite 
wave-number, i.e., the conductivity of an anisotropic current. 

The outline of this paper is as follows: In section 2.2, we describe how to introduce 
the background quantities in the defect system and describe several problems that arise in 
the D3-D7 setup - which motivate us not to pursue the massive D7 case. We then show 
in section 3, how to obtain the conductivity from linear response theory in our case, and 
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also demonstrate the effects of electromagnetic duality at finite density and finite magnetic 
fields. In section 4, we then derive analytic results in various limits, first in the isotropic 
DC limit u> —> and in the small frequency expansion beyond the DC limit. Then, we 
consider the small temperature limit q » T at small frequencies u <C q in various regimes 
of the density and magnetic field and finally we obtain the diffusion constant and "electric" 
permittivity and consider the hydrodynamic regime (u, q <C 1). The numerically-obtained 
full spectral curves are presented in section 5, where we also later present and discuss 
the purely dissipative poles that we obtain numerically and finally study the quasiparticle 
poles in the correlator that we obtain numerically both directly and from the spectral 
curves. We present the explicit form of the induced metric on the brane in appendix A, 
and in appendix B, we review some basic properties of weakly coupled systems in order 
to introduce the terminology and remind the reader of the reader of some intuition and 
generic expectations. 



2. Turning on the Background Parameters 
2.1 D3 = (N=4 SYM) background 

Let us remind ourselves of the super Yang-Mills background. We start off with the well- 
known AdS§ x S 5 background of N c D3 branes in the decoupling limit corresponding to 
an N = 4 SYM theory on the boundary [1, 2, 3] with U(N C ) gauge group. We work in the 
limit of N c — ► oo, at Yang-Mills coupling g\ M = 2irg s — ► in the field theory, such that we 
consider the large t'Hoft coupling limit A = gy M N c — ► oo and we can use the supergravity 
limit as L 4 = 4irg s NJg — > oo. At finite temperature T = the background metric is 
written as 

ds 2 = £ (-h(r)dt* + dxl) + ^ + r 2 ^|) , = . (2 . 1} 

Considering only T > allows us to go to dimensionless coordinates u = ^, t = x = 



ds 2 = ^ (-(1 - u*)dP + dx% + + u 2 d^ . (2.2) 

In the field theory side, all fields in the SYM theory transform in the adjoint repre- 
sentation of the gauge group. The most well-known approach to introducing matter fields 
transforming in the fundamental representation is to insert probe D7-branes in the radial 
direction into the supergravity background [5], intersecting in the "flat" (SYM) directions 
with the D3-branes. If the intersection overs only part of the flat directions, this will cre- 
ate a defect field theory, where the fundamental fields are only supported on a subspace 
within the four-dimensional spacetime of the gauge theory. As in [28], we will consider 
a (2 + l)-dimensional defect by inserting Aj Dp-branes, with three dimensions parallel to 
the SYM directions and p — 3 directions wrapped on the S" 5 , considering both probe D5- 
and D 7- branes. If we consider the supergravity background as the throat geometry of A^ c 
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D3-branes, our defect constructions are described by the following array: 



0123456789 



background : D3 x x x x 

probe : D5 x x x x x x 

D7 xxx xxxxx 



(2.3) 



The D5-brane construction is supersymmetric and the dual field theory is now the 
SYM gauge theory coupled to N f fundamental hypermultiplets, which are confined to a 
(2+l)-dimensional defect. Note that the supersymmetry has been reduced from J\f = 4 to 
N = 2 by the introduction of the defect. In the D7-brane case, we have lost supersymme- 
try altogether and the defect supports N t flavours of fermions, again in the fundamental 
representation [24]. One should worry that the lack of supersymmetry in the latter case 
will manifest itself with the appearance of instabilities. In [28], we showed how this can be 
avoided, an we discuss in section 2.4 how this instability becomes apparent in the scaling 
dimension of the scalar field that corresponds to the deformation of the S 4 of the D7- 
worldvolume inside the S 5 background. There, we allso discuss some problems related to 
the reliability of the quenced approximation that we consider in this paper. In this limit, 
N f <C N c , the D5-branes may be treated as probes in the supergravity background, i.e., 
we may ignore the gravitational back-reaction of the branes. For the D7-branes, however, 
this is only true locally and not in the asymptotic regime. 

As we commented above, a similar holographic framework has been used extensively 
to study the properties of the N = 2 gauge theory constructed with parallel D7- and 
D3-branes, i.e., the fundamental fields propagate in the full four-dimensional spacetime 
- e.g., see [7, 6, 8]. There, it was found that if a "quark" mass M q is introduced for 
the hypermultiplets, it was found that the scale M fun ~ M q j\f\ plays a special role in 
this theory. First, the "mesons", bound states of a fundamental and an anti-fundamental 
field, are deeply bound with their spectrum of masses characterized by M fun [9] . Next at 
a temperature T ~ M fun , the system undergoes a phase transition characterized by the 
dissociation of the mesonic bound states [8]. The meson spectrum is characterized by 
the same mass scale M fun [29, 30] and these states are completely dissociated in a phase 
transition at T ~ M fun . A similar bhaviour can be observerd for defect in the D5 case, that 
is T-dual to the mentioned D3-D7 configuration. We will discuss this briefly in section 2.3 
and more in detail in [31], where we also find some properties that are new in the case of 
a defect. 

2.2 Introducing the defect 

In the supergravity limit, the D5 brane action of the U{\) subgroup of the U{Nf) is just 
the DBI action plus a Chern-Simons term 



where the factors of Nf arise from taking the trace over the flavor degrees of freedom, 
arising from the stack of Nf coincident branes. To simplify things further, we work in the 



S = -T 5 N f [ y/-det(P[G] + 2irP s F) + T 5 N f [ C (4) A 2nl 2 s F , 



(2.4) 
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quenched approximation Nf <C N c , such that we can ignore the backreaction of the probe 
branes. 

Preserving translational invariance in the fiat directions and rotational invariance on 
the sphere, together with the choice of the embedding (2.3) dictates the induced metric on 
the D5 brane to be of the form 

ds 2 = g (-(1 - u*)dP + <& 2 + (l + (1 - u 4 ) (z> { u) 2 + u 2 ^^) ) ^ 

+ u 2 {l-^{u) 2 )dn 2 2 ) . (2.5) 

In a previous paper [28], we discussed the case of the trivial solution F = for the 
U(l) background in the fiat directions and for the S 2 radius * = 0. Even though 
has a tachyonic mode, corresponding to shrinking to zero size, it's mass lies above the 
Breitenlohner-Friedmann bound my = 2 2 > — [34], such that it does not cause 

AdSi L 4AdS 4 

an instability. The discussion of the setup was limited to the case of vanishing "quark" 
mass for the matter on the defect, vanishing net density of matter and antimatter on the 
defect (net baryon number density) and vanishing external magnetic field applied to the 
defect, but a flux F = j^df^ on the compact sphere was turned on. This corresponds on 
the gravity side to having an extra set of q D3 branes pulling on the D5 from one side of 
the defect and on the field theory side to having an extra number of colors, 5N C = q on 
that side of the defect. Both from the embedding geometry z(u), and from the resulting 
quasiparticle spectrum in the field theory, it was argued that this flux also introduced a 
finite width, Az, of the defect. The embedding was found to be 

= 7TW (2 ' 6) 

which has in principle some analytical solution. Here, as everywhere in this thesis, we use 
the notation (■)' := d u (-). 

2.3 AdS/CMT Dictionary (Supersymmetric Case) 

Now the situation is slightly more non-trivial, as we wish to introduce finite values for the 
mass, baryon density and magnetic field. Using the AdS/CFT dictionary in [3] in analogy 
with the 3 + 1 dimensional system, e.g. [26, 49, 11], we find the gravity dual of the baryon 
density 

p = -2tt 2 -— = 2TT 2 (TTT) 2 e lim A'iu) =: -27r 2 (7rT) 2 lim E(u) (2.7) 

and magnetic field B to be related to a non-trivial U{\) background on the brane: 

F| u ^ = -E(u)dt Adr + Bdx A dy =: F E + F B . (2.8) 

We can also define the (asymmetric) background metric 

G = g + F. (2.9) 

By analogy with the 3 + 1 dimensional D3 — D7 system, we can repeat the arguments in 
[8, 26, 27, 7], and associate a non-trivial embedding ty(u) with a finite quark mass M q and 
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dual condensate C. This condensate has on the one hand an interpretation as a chemical 
potential for M q and on the other hand is considered in QCD contexts considered as the 
order parameter of chiral symmetry breaking. 

We can parametrize the S b as dfi 2 = dip 2 + cos 2 tp dVL^ + sin 2 tp dfl 2 . and put the D5 
on the first S" 2 , such that sin^ =: \£ gives us a suitable scalar describing the size of the 
sphere. In the parametrization (2.8) the DBI-CS action becomes then 

S = 4nL 2 T 5 J da 4 (V-detGVC 1 " ^ 2 ) 2 + P + fu A d u z^j (2.10) 
= 4ttL 2 T 5 J da A (V^det^^l + F 2 ^1 + F| V(l - * 2 ) 2 + P + fu A d u z^j , 
and one trivially finds the background solution 



E{u) 



d u z 



B = const. (2.11) 



^/l-*(u) 2 ^l + (f 2 + (p 2 + B 2 ){1 + f 2 ))u* + (1 + B 2 u^{u) 2 (^{u) 2 - 2) 

-fy/l - 2 (it) + u 2 h{u)^\u) 2 
y/l - ^{u) 2 ^l + (p + (p 2 + B 2 )(l + p))u A + (1 + B 2 u 4 )y(u) 2 (V(u) 2 - 2) 



for all the physically relevant fields, except for *(«), because that one enters the action 
both directly and with one derivative. As a reminder, our notation is h(u) = 1 — u 4 . 
For convenience, we defined the dimensionless parameters p := 27 r 2 (^T) 2 £ and = (ttT) 2 • 
We see that the width of the defect from the brane picture, z max := lim u _>i z(u), de- 
creases as we increase B and po as they appear only in the denominator. This may appear 
somewhat counter-intuitive from a weakly coupled point of view, but it is what we should 
expect, as the system is strongly coupled, or the correlation length diverges, and hence the 
"contractive force" scales with the total number of particles. 
The equation of motion for 'I'(u) becomes 

2(1 + BV)(1 - ^ 2 ) 3 + u 2 (l - u A ) (1 - (f + (p 2 + B 2 )(l + f))u A + (1 + £ 2 u 4 )* 2 (* 2 - 2))^ /2 



4 (1-^ 2 )^(1-^ 2 )(1-^ 2 + {u 2 - u e )V 2 ) (1 + (/ 2 + (p 2 + B 2 )(l + / 2 ))u 4 + (1 + B 2 u^ 2 ^ 2 - 2)) 

(2.12) 



. , 1 - n 4 l + (P + (p 2 + B 2 )(l + P))u 4 + (1 + B 2 u 4 )^ 2 {^ 2 - 2) 



u 2 V (l-^ 2 )(l-^ 2 + (u 2 -u 6 )^' 2 ) 

which has no analytical solution, except for some limiting cases. For u — > 0, it is easy to 
see that the solution becomes 

* ~ fnu + cu 2 , (2.13) 

where rh and c are dimensionless free parameters that are determined by the boundary 
conditions. Now, we see that the argument of the T-dual case of the D3-D7 system with 
3 + 1 intersecting directions [8, 26, 27, 7] also applies to our case, and the quark mass M q 
and condensate C are given by 

M i = SPr? = and C = V24ir 2 cr 2 N f l 2 T 5 = cT 2 N f N c . (2.14) 

Z^'^TTlg Z" 5 ^ 47T 
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This can be straightforwardly obtained from the results in [8], but also in our case we see 
that this relates to the length of a string spanning on the sphere from the D3 branes to 
the D5 branes. In [31], we discuss this more in detail and verify that C is indeed the dual 
chemical potential to the mass. 

In order to find the solution for the full geometry for a given mass however, we need 
consider the equation near the horizon, where (2.12) reduces to first order, 

*'i , = 1 (i + 5 2 )tto(i-tt8) 2 (2151 

' 2(1-M/2)2 + /2 + /5 2 + j B 2 ( 1 + (1 _^2 )2 ) ' ^ > 

effectively relating rh and c. The only remaining boundary condition at the horizon is 
then ^ u ^i = ^o- Because this boundary condition is dictated to us by the equation of 
motion, which becomes first order at this point, and so is also dictated to us the other 
boundary condition \l/|u=o = 0, we use rh instead as a boundary condition, but we have 
to find recursively fo r a given value of rh. This is because implicitly, on-shell, c is a 
function of rh and starting to integrate at u = with some random combination of c and 
rh means the equations of motion cannot be on-shell as we approach the horizon. 

At vanishing density po and vanishing compact flux /, we find again that the black 
hole embedding which gives us free quarks, limits M q < Merit- At M cr u, we have the 2 + 1 
analogue of the phase transition that was found for the 3 + 1 system in [8, 7]. It turns 
out that the critical mass decreases as we turn on the magnetic field in the 2 + 1 field 
theory. This is discussed in detail in [31]. This phase transition disappears (at least in 
our case where we consider only the U(l) background) as we turn on either a finite baryon 
density po, or as we choose the compact magnetic flux / to be non-zero. Essentially, this 
happens because the charge can only be supported by the blackhole embedding and the 
action becomes singular when *(u) = at finite flux just as in the 3 + 1 dimensional case 
in [8]. 

Some of our studies of the effects of finite masses, we will have to consider finite po or 
/, to allow for sufficiently large masses. In the limit of very large masses (i.e. — > 1_ near 
the horizon), one can see that over u g]0, 1], the equation of motion for \l/ is also solved by 
\l/ ~ 1 - to see this in eq. (2.12), one needs to take ^f' — > first. This demonstrates how a 
new length scale arises for large masses, m > 1 as the profile splits approximately into two 
parts - one with u > u m ~ i and ^ ~ 1 for some value u m and one given approximately 
by the asymptotic solution (2.13). Around u m , V diverges. It would be interesting to see 
whether this limit has any relation to the recent discussion of non-relativistic AdS/CFT 
[16]. 

A more thorough discussion of the thermodynamics and the phase structure can be 
found in [31]. 

2.4 Non-supersymmetric D3 — D7 intersection 

The non-supersymmetric case is very similar to the D5 case as it differs in the massless 
case only by the geometry and field configuration in the S 5 factor. Now we parametrize 
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the S 5 in the bulk space as dO 2 = dtp 2 + cos 2 tp d£l\, such that we have the induced metric 

ds 2 = g (-(1 - „*)dP + dil + (l + (1 - u 4 ) (z>{u? + ^^p) ) ^ 

+ u 2 (l - ^(w) 2 )c«7 2 ) . (2.16) 

and we set up an instanton on the S 4 instead of the magnetic charge on the S 2 . The 
coupling to the five-form flux comes now via the Chern-Simons term 

^f^T 7 N f f C« A F A F . (2.17) 

2 JD7 

This CS term however also causes the D7 setup to differ from the D5 setup in the massive 
case, as can be seen most easily by integrating this term by parts to give us (modulo a 
total derivative) 

(27r ^ )2 T 7 iV f f F^AAAF 8T 7 iV f ^ du , = : / iAF. (2.18) 

The term on the right hand side with the factor \P' arises from the fact that the deformed 
embedding of the S 4 inside the S 5 causes the dual F 5 on the sphere to pull back to the 
brane worldvolume. Integrating out the S A of the worldvolume gives us then the right hand 
side, which is just a Chern-Simons term with radius-dependent coupling. This was used 
in [32] to obtain a Hall effect in a setup that would be considered from the perspective of 
this paper unstable. This term will obviously modify the two-point functions, and will be 
interesting to consider in further work, but there are some problems with the massive D7 
case that we will outline below, so we will here only consider the massless D7 defect. 
The instanton solution was found in [33, 35] and is outlined also in [28] and yields 

I TrF A F =: q 7 G Z (2.19) 
8tt z J S 4 

and the 5 4 factor 

I d 4 n = ^ (N f L\\ - * 2 ) + 6vr 2 ^ 4 |g 7 |) , (2.20) 
Js 4 & 

such that the Ansatz (2.8) puts the action into the form 

^-N f L A T 7 j da A (7detG((l - * 2 ) + \Q\) + Qu A d u z^ (2.21) 



S = 



where we defined Q = 67r 2 |%|^; = ^f-J^- Having an isotropic solition on the S A , i.e. 

1-1 1 II. N f( N ^~ 1 ) T 1 1 

having the corresponding symmetries unbroken requires g 1 > q 7 . In the general 

case (\P 7^ 0), this action gives different solutions than (2.11), however in the case \& = 
the solutions are given precisely by (2.11), provided we replace the flux parameter with 



Q_ 



In [28] it was found that the mass of the tachyonic mode of the S radius of the D7 
probe brane satisfies the BF bound only for f 2 > 49/32, and a quick calculation shows 
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that this also happens to apply in this background - independent of p and B. This, we will 
see, is reflected in the asymptotic behavior of ^(u). We will not bother the reader with 
the lengthy form of the equations of motion for ^, however, we note that the asymptotic 
solution takes the form 




*(» - »" = ■ ni - ± t 2 + 4Q (2.22) 

which implies that above the BF bound, the solution will be a power law, and below the 
BF bound it will be oscillatory - indicating the instability. Above the BF bound, we could, 
in principle, identify those two modes with a "mass-like" operator (and a "condensate" 
operator) of non-integer conformal dimension, motivated by the fact that this is related to 
the separation of the D3 and D7 branes in the sphere. Possibly one could interpret this 
behavior with a "running" mass. However, there is a very significant problem that may be 
more worrying than the instability at vanishing mass: A solution of Nj D7 branes causes 
an asymptotic deficit angle of ^ from backreaction (see e.g. [36]). At the necessary finite 
fi, and hence q 7 of order A, we need at stack of D7 branes as determined by the limit 
Nf ^ N Q f - > q 7 for spherically symmetric solution from [35]. In order to still have the same 
kind of field theory with the same symmetries, we require Nf 3> 12, which implies that the 
solution cannot connect to an asymptotic space-time. Note that this was exactly the S 4 
factor which broke supersymmetry, and this factor would be highly modified for solutions 
with q 7 below this bound. Below the BF bound the oscillatory asymptotic solution is 
non-physical in our setup and furthermore the corresponding operator would have complex 
conformal dimension. It seems that there are some non-standard ways to interpret the 
ff = case in the context of the quantum Hall effect [25], but we will not pursue the 
^ 7^ case in the rest of this thesis. It is a noteworthy curiosity, that in the absence of the 
pullback of the CS term to the flat directions, the resulting spectral functions are identical 
to the D5 case under an appropriate identification of the mass-like operators. 

3. Computing the Conductivity 

One common approach to compute the conductivity that we will pursue here is linear 
response theory, i.e. applying the Kubo formula 

ir~ l T(Jij = &ij = —Cij or for convenience &ij := irTaij = —Cij , (3.1) 

LO LO 

that gives the conductivity for currents resulting from small perturbations in terms of the 
retarded Green's function, which is given in terms of the correlator as 

C ij (x-y) = -i0(x°-y o )([J i (x),J j (y)]) . (3.2) 

We define for later convenience a = ^ . Since the baryon number current Jj is dual to the 
gauge field of the U(l) subgroup of the U (Nf), Ai, the correlator is given by the variation 
of the on-shell action 

5 2 S 
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where Aq is the boundary value of the gauge field at the asymptotic boundary u = 0. The 
action for the gauge field becomes to second order 



(3.4) 

where G is the asymmetric combined metric (2.9), and g\ is defined as 

1 = 4vr(2^ 2 ) 2 L 2 7lTrT 5 = V^Tp-^= • (3.5) 
<?4 K V 2A 

\/(l-^ 2 ) 2 +f 2 

In some sense, there is now a radius dependent coupling — — ^/ 1+ p — ' a ^ wa y s g° es to 
unity asymptotically or obviously everywhere in the massless case. Surprisingly, the gauge 
field background dies off sufficiently fast asymptotically, such that as u — ► 0, the action just 
becomes the Maxwell action with coupling and in a suitable gauge A u = the correlator 
is still given by the asymptotic mode function 



C ij = e 5 (d u Aj) 



VT+P SAi 



(3.6) 



7rT 



where eo is defined as eo = 

94 

3.1 Electromagnetic duality 

In this background, we see that the effective action for the gauge field is not invariant 

e 2 u 2 

under electromagnetic duality F — > *F. Hence, the relation C xx = — -A — that was found 
in [10] does not apply in this case. Since the DBI action at constant coupling, i.e. in the 
massless case, however still obeys this duality, one would expect that it survives in some 
form under the exchange of the magnetic and electric charges on the probe brane, i.e. under 
the exchange of the density and magnetic field in the field theory side. To quantify this 
further, let us look at the transformations under F — > *F of the Fourier-transformed gauge 
field in the gauge A u = that led to (3.6), obviously at M q = 0. The relevant components 
are at asymptotic infinity: 



(*F) tx \u=o = -V 1 + f 2 F uy \u=o and {*F) ty \ u=0 = + f 2 F ux \ u=0 , (3.7) 

such that the variation w.r.t. the gauge field becomes in terms of the transformed gauge 
field, denoted in abusive notation as *A: 



d 




(9A' y d 


dA x 
d 




\8A X d(*A) x 
(dA' y d 




dAy 


u=0 


\8A y d{*A) x 



dAL d \ d /„ d _ 1 



dA X x d(*A) y ) u -eoV xy d(*A) x ~° xx d{*A)J u=Q (3 ' 8) 

OA' d \ d ( 1 d \ 

dt y d^A)y) u = V, {^d^AT x + ^d^ATj u=0 ■ (3 ' 9) 
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Rewriting the conductivity obtained from (3.3) then in terms of the transformed fields 
gives us 

&xx = ~2 xy@ xy\*<j\xx + d- xx a xx [ka}yy) (3.10) 
e o 

**yv = T2 (^yy^yyl*^]xx + ^xy^xy[*S]yy) (3-11) 
e o 

&xy — ~~ 2 i&xy@yy[* : @\xx ~\~ & xx &xy\*G\yy ~\~ (&xy&xy @xx®yy) \*0\xy) (^ - 1^) 

where we used (3.6) and S[F] = S[*F], and defined [-ka]^ = s[*A] S js[kA]- - Since F <-> *F 
exchanges the electric and magnetic charges on the probe brane defined at infinity - i.e. 
exchanges density and magnetic field in the field theory side - *cr is just the conductivity 
and the exchange of p and B. Finally, we can solve for *cr and obtain 

where a, b G {x, y}. 

This result is remarkable, since it relates the transport properties under the exchange 
of two quantities that are completely distinct in nature from the condensed matter point 
of view. Furthermore, it applies to a whole class of strongly coupled 2 + 1 dimensional 
systems, whose gravity dual obeys the electromagnetic duality. Hence, such a relation 
is a generic prediction from AdS/CFT for a quantum critical 2-dimensional system. For 
theories not obeying (3.6), there may potentially be additional terms in (3.13). It seems 
that this is an implication of the "particle- vortex duality" found in [43, 44], extended to 
finite frequencies and accordingly to a complex conductivity tensor. Certainly, this duality 
does not generate the full SX(2,Z). A candidate for the second generator is simply a shift 
in the theta angle and the corresponding Hall conductivity from the Wess-Zumino term 
in the action that was outlined in [28]. For another discussion of 5*1/(2, Z) symmetry in 
the context of AdS / CFT in a somewhat different limit that appeared very recently on the 
arXiv see [45]. 

This duality holds always in the massless case to numeric accuracy. Hence it is not 
possible to visually "compare" the result in a plot. 

3.2 Explicit Computations 

To proceed further let us start by writing out the equations of motion explicitly: 

A y Q = {\f^GG yy G uu A' y )' + s[^GG vv {G xx dl + G tt d?)A y 

+ (^GG tu G xy )' (d x A t - d t A x ) (3.14) 

A u a = G tt d t A' t + G xx d x A' x (3.15) 

A t Q= (V^GG a G uu A' t )' + V^GG tt G xx (dlA t -dtd x A x ) (y^GG tu G xy ^j d x A y (3.16) 

A X Q = (y/^GG xx G uu A' x )' + V^GG U G XX (d?A x - d t d x A t ) + [^GG tu G xy ^j d t A y .(3.17) 
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For convenience of the reader, we stick here to the concise notation in terms of G and 
summarize the exact form of the components G in appendix A. Also, in this expression, and 
for the rest of this thesis, we absorbed the radius-dependent coupling into the determinant 
of the metric, in somewhat abusive notation: 

V{ ' W-det G V^G . (3.18) 

Finally, we also remind ourselves that G yy = G xx , so while keeping them distinctively 
for didactic reasons in most places, in some places they will be interchanged to simplify 
expressions. 

It can be easily verified by using the equation for A u , that the equations for A t and A x 
are degenerate. Hence, our strategy will be to eliminate A t from the equation for A x and 
produce an equation for A' x , by multiplying (3.16) with y ' —GG yy G uu and differentiating 
with respect to u. This gives us: 

0= (y^GG yy G uu A' y y + ^V^GG yy G xx dl + V^GG yy G u d^A y (3.19) 

{{^GG^G xy )') 2 (y/^GG tu G xy )' , 

+ V^GG U G XX y + \[^GG tt G xx x 
0= (y^GG yy G uu A^ + (y^GG yy G xx d 2 x + y/^GG™ G* 6%) A x (3.20) 

- (-GG«G-G~G«y A , + ( _ a)GttG ^ G „((-^GG'«G'*)' V 

^QQttQxx x \ ) y ^GG tt G xx y J 

where A x = — . These equations separate at vanishing density or vanishing 

magnetic field - as they should, because we do not expect a Hall effect in this case. 

Using (3.17), we can recover Ar| M= o = —i ^Q-vp) ; which will allow us to compute the 
conductivity tensor. Near the horizon, the solutions become approximately 

A y = ( i- u r/ 4 (4°)+4 i )(i-n)) , a x = (i-ur/ 4 (4°)+4 i )(i-«)) , (3.21) 

where A y ^ and A x ^ are arbitrary constants and Ay 1 " 1 and are determined straightfor- 
wardly in terms of pQ, B, ^q, A y °^ and A x °\ but somewhat lengthy and without physical 
insight and hence omitted here. To compute the conductivity, we could then fix the bound- 
ary conditions for A y £ { 1 , 0} and A x G { 1 , 0} at u — ► and enforce the leading behavior 
of (3.21) as a boundary condition at u — > 1 as done in [28]. This is however a numerically 
non-trivial boundary value problem. Hence, it is more reliable and less time-intensive to 
simply enforce (3.21) for two independent choices of {A y °\ A^}, labeled a and b, to ob- 
tain {Ay(a),A'y(a),A x (a),A' x (a)} and {A y (b), A' y (b),A x (b),A' x (b)} at u ->■ 1 and then use 
the linearity of the problem to compute the variation in (3.6) exactly. Furthermore, this 
strategy is very suitable from a computational point of view, as it allows us to naturally 
the parallelize solving the equations of motion, i.e. the most time consuming step, on a 
dual-core processor. 
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Putting everything together, we finally obtain: 



a = s 



2- A x {a)A y {b) 

QJ 



-A x {b)A y (a) 



/l | p A y (h)A y \ a )-A y (a)A y '{b) 



1 



A x (b)A x '(a)~A x (a)A x '(b) 



A x '(a)A y (b)-A x '(b)A y (a) 

iyjb) Ay' (CI) -Ayja) Ay' (b) 

A x '(a)A y (b)-A x '(b)A y (a) ^pu, A x '(a)A y (b)-A x '(b)A y (a) 



A x '(a)A y (b)-A x '(b)A y (a) 
A x '(a)A y '(b)-A x '(b)A y '(a) 



Formally, this is asymmetric, such that the (anti)symmetry of the numerical result is a 
check for the consistency and the accuracy of the numerical solutions for A x and A y . We 
also note that in the limit of {po, B, = 0, we just recover the equations that were found 
in the conformal case in [28]. In principle, the duality from 3.1 suggests that there may 
exist a field redefinition for A y and A x , such that the asymptotic solutions for the resulting 
fields are exchanged under p <-> B. However, there is no guarantee that this redefinition 
can be written analytically. 



4. Analytic results 

4.1 Isotropic perturbations with small frequencies 
4.1.1 DC Limit 

One limit of obvious interest is the isotropic DC limit (i.e. ui, q <C {1, p~ }). To 

analyze this case, we define a new radial coordinate s, g£ = (\^GG yy G uu ) and study 
the original equations of motion. Now, they just read 

A y : = d 2 s A y + V^GG yy G uu (V^GG yy G xx d 2 + V^GG yy G tt d 2 )A y 

+ s/^GG yy G uu (y^GG tu G xy Y (d x A t - d t A x ) (4.1) 
A u :0 = G u d t d s A t + G x d x d s A x (4.2) 

A t :0=^d 2 A t + ^GG yy G uu [^j d s A t + (-G)G yy G uu G u G xx (d 2 A t - d t d x A x ) 

- V^GG yy G uu [^GG tu G xy ^j d x A y (4.3) 
A x :0= d 2 s A x + (-G) G yy G uu G tt G xx (d 2 A x - d t d x A t ) 

+ V^GG yy G uu (y^GG tu G xy ^j ' d t A y . (4.4) 

First, we consider the equations in the limit u — > 1. In this limit, we have §r = , . i+b 2 +f 2 == 

du ^(l+/ 2 )B 2 (/ 2 +(l-*g) 2 )+(l+/ 2 ) 2 ( 

0(1) =: + 0(1), such that s ~ -c s ln(l - u 4 ) or (1 - u 4 ) ~ . Note that in this 

coordinate, the horizon is located at s — > oo and asymptotic infinity is at s = 0. Now, the 
equations of motions reduce to 

o - + (1 + ^ ^ + * - gg^±ig + /2 + " - (4,, 

= a A t , = d 2 A t , (4.6) 
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up to order e c « . This has the solution 

A t = A° t , A x , y = Al y e™ ; (4.7) 

Jb 2 (p + (i - + (i + p){~ P 2 + f 2 + (i - ^) 2 ) 

f = iiv 1 + / 2 - 



1 + B 2 + / 2 

Then, we consider the region of s £ [0,0(1) x c s ]. To obtain the diagonal and Hall 
conductivities, we set e.g. A x = and At = at s = 0, and study the linear depen- 
dence of d s A y and d s A x on A y . Combining the equations (4.1) to (4.4), we find that 
d 2 A y < O ((a; 2 + q 2 )(l + f 2 + B 2 + p 2 )j A y , such that the change in d s A y over this region 

,2 , ~2\„ n i j-2 i d2 i ~2\ 



is 5(d s A y ) < O [{id 2, + q )c s (l + f 2 + B 2 + p 2 )j A y . Hence, we have to leading order 

d s A y = ivA y , (4.8) 

which gives us the isotropic DC diagonal conductivity and also allows us to obtain the Hall 
conductivity. 

Looking at A x , we find that to leading order in Co,q, we have 



= d u d s A x + iuAydu [^GG tu G xy j or (4.9) 
= d 2 A x + iuA y V^GG uu G yy d u [^GG tu G xy ^j , 

and hence A x is a consistent solution near the horizon at large s. In the asymptotic region 
at small s and in the limit of Co <C 1 the first integral of u can be done analytically, such 
that we obtain 



d s A x = -iCoA° y 



-GG tu G xy 



(4.10) 



J u=l 

since we should have by consistency d s A x — > as s — > 0. Hence A x \ u= o ^ A y x 0{Co) <C 

^4.j/|u=0 

Finally, we note that = 1 as it — > 0, such that we can write down the DC conduc- 
tivity tensor 



a —eq 



lB%p+(l^lY) + (l+f 2 ){P 2 +P+(l-^l) 2 ) Bp 



l+B 2 +f 2 l+B 2 +f 2 

Bp ^(/2 + (l-<l/g)2) + (l+/2)(p2 + /2 + (l-^)2) 

1+B 2 +/ 2 1+B 2 +/ 2 



It is straightforward to verify that at \&o = 0, i.e. M q = 0, this expression obeys the duality 
under exchange of a <-> ^o' -1 , x ^ y, p ^ B. Obviously, to obtain the full dependence 
at finite mass, we have to invert M q (f,p,B,^o), in order to obtain \&o as a function of 
{M q , f, p, B} - but this is not possible in closed-form. At small masses m <C 1 and small 
/, p and -B, however, one can use rh ~ $?o, and at large quark mass m > 1 at finite p and 
/, we have ~ 1 and the result becomes independent of the quark mass. 

If we compare this result with the Drude conductivity (B.2), we find that this is 
qualitatively what one would expect. We can identify ^J^~p = UJcT, > z"^ 1 = IfrP an< ^ 
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/in = jqTp . The tact that the expression tor fin is 

somewhat complicated is not surprising, since it results from the density of quark-antiquark 
pairs in thermal equilibrium. What is somewhat surprising is the fact that at finite /, there 
is only very limited dependence on the quark mass M q - because one might have thought 
that (at vanishing po) n is strongly suppressed at large M q - but one should not interpret 
too much into this result. What comes as expected though is the fact that //n oc p at large 

4.1.2 Small frequencies 

Next, let us try to extract the subleading terms in the conductivity at small frequencies. 
To do so, we perturb the equations of motion for A y and A x (4.1), (4.4) around the DC 
solution by taking A X;V — ► A xy e lus + A$ y . The equations of motion for A$ y becomes then 
at q = 

& a Ag) = 5{Q tu G uy ) u 2 A^e ivs + iojg uv (g 1 ^' A y e ivs , (4.11) 

s%aW = s(g tu g uy ) t?Afe ivs - m uy (g 1 -^ aJ vs , (4.12) 

guy ._ ^GG yy G uu , G ty := V^GG yv G tt , G 1 := V^GG xy G tu , 
S(-) := (•) - (■)«-»! • 

We also write out the symbols g in appendix A. For simplicity, we choose as above 
s\ u =o = 0, and we use u as a variable to work with. Before proceeding, we look at the 
correction to the conductivity: 



°™ ~ ~ £o uSA y 



-DC I 1 1 1 °s^y _ siy \ 
yy { + iv a (0) ' A {0) / «» 



Primarily, we are interested in the 0(lo 2 ) corrections to the real part of the conductivity, 
so we need to keep track of A^ up to 0(lu 2 ) and d u A^ up to 0(lo 3 ), which coincides with 
the accuracy of the first perturbation, as the natural expansion parameter is lo 2 . In the 
case of pqB = 0, only the diagonal term in the equations of motion contributes, so we find 
to the relevant order 

du(l + iM«)) g„, 1 (4.13) 
jV =S ?W]l*±jl*Wgp-. (4.14) 

We remind ourselves that s(u) = J™ g^, such that s(u) and the first integral of ^ g uy — - 
can be easily computed analytically at M q = and expressed in terms of hypergeometric 
functions. The second integrals however have to be computed numerically even in the 
massless case. To demonstrate convergence, we note that any combination of the form 
^GQHv G af3 ig finite a t u ^ o and at it ^ 1 we have s oc - ln(l - u 4 ), g uy oc (1 - u A ) oc 
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5 (Q tu Q uy Y The convergence of J duln(l — u A ) is also the reason why we could expand the 
exponential at sufficiently small to. 

Including the case of B x po / is slightly more tedious. First, we compute A x up to 
Oiio 1 ) under the condition that A x \ u= q = 0. To do so, we first need to integrate d u (d s A x ) 
in (4.4). The condition Aj.| M=0 = implies then that d s A x \ u= \ ~ iv J^~^ d s A x , such that 
we can, in the limit of small w, use d s A x \ u= i = 0. Hence, we get: 



A, 



iuAf 



guy 



y Jo V ^ uy $ uy Jo $ uy 



vuAf 



(4.15) 



io <3 uy 

Here, as in the rest of this section, we performed the integration by parts in order to limit 
the number of consecutive integrals to two integrals. Now, we can compute the additional 
contribution to Ay~\ which can still be written in terms of double integrals, with the first 
one computable analytically at M q = 0: 



£ 2 A (0) 



guy 



I 

J u 



du- 



guy 



(4.16) 



5AM 



iu 2 uA^ 



g : - / du \ s (u) 



guy 



I 

J u 



dug 1 ' s{u) 



guy 



guy 





l 


guy J 


guy 


— — / du 




Q uy Jo 


guy 



guy 



cD 2 A (0) 



guy 



J u 



du 



guy 



(4.17) 



The integral for SAy 1 " 1 might seem divergent to the reader, but by close inspection it is 
apparent that the integrand is finite as u — > 1. Finally, we can write the correction to 
the diagonal conductivity which simplifies significantly after some simple algebra: After 
setting u = 0, we can eliminate the first term in the first line and all of the second line in 
(4.16) and then it turns out that most of the terms in A y 1 ^ and d s Ay^ are pairwise equal, 
such that we obtain 



a, 



yy 



f DC 

yy 



- 2 / 

J 



du 



6{G 



1 sigtygvy) r l 5(G L ) 2 

dus(u) K niii ' - 2 / dus(u) {y ' 

Jo 



guy 



guy 



guy J Q guy 



( r 1 5(G ty G uy ) f 1 
+ ieo& / du + / du 

[Jo G u y Jo 



duG s(u) 
guy I 



guy 



(4.18) 



For completeness, we can also compute the contribution to the Hall conductivity. To 
do so, we again consider a pertubation that keeps A| u= o = 0. The Hall conductivity will 
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then be to order lo 2 



DC 

xy 



A 



(1) 



A 



(o) 



+ 



(1) 



0.4 

(o) 



(0) 



(4.19) 



u^0/ 



We already know — ^y, so we only need to compute d s A y x ' . There will be two contributions, 

Ay 

from the diagonal and off-diagonal terms in the equation of motion for A x . Using as zeroth 
order the first term 

(4.20) 



iuA^ 



du- 



d,AW 



(4.21) 



/o G uy 

from (4.15), we find that the contributions from the diagonal term in (4.11) is 

guy J Q guy y J Q guy 

where the second term comes from the oscillatory behavior towards the horizon at large s. 
In the off-diagonal term, let us first write the O(oo 2 ) term 

"guy J 

ru du 



Sd s A^ = CjvA^ (g^s{u) - J du- 



8A<£> = uvAW , 
lo Q uy 



1 G ± 
du— — 

du G uy 



(4.22) 
(4.23) 



giving rise to an 0(<jj 3 ) term 



<5i&4« = 2iCou 2 A\^ Cdu 

Jo 



g ± s(u) 



(4.24) 



io G uy 

Using the fact that ^- l | u= q = 0, the direct 0(oj 3 ) contributions from the cross-term 
are read off from (4.13) and (4.17): 



S 2 d s A^ = iuis 2 A^g ± s(u) 2 + 2 / du—{u)Q h 

Ju s 

/ guy I- 



iu 3 A^ g 



\_5{g tx g uy ) 

guy 



+ 



Cdu^- [ 
Ju U Q UV Ju 



i 5{Q tx G u y) 

du 

guy 



_ e x r* i(5 x, r^m: 

J„ e«» 1/ ' y G"« y 
J„ G'"< \ Jo e*» A 



guy 



(4.25) 



j guy 

We can note that the first term in each line vanishes if we take u — > 0. In this case, 
also the contribution from (4.17) and the last two "sub-terms" all combine into one term, 
the contributions from (4.24) and the first line are equal, as are the second line and the 
contribution from (4.21). Hence, we the result can be written as : 

pB 



'xy 



1 + f 2 + B 2 



+ uj £o / du 



Jo 



f 1 5 i<3 

2iveQ I du 
Jo 

5(g tu g x y) - 5{g^ 



guy 



guy 

Jo 6"» ' Jo S" Jo 



du 



(4.26) 
guy 
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Let us now look at the fruits of this algebra. 




ciio- 



Figure 2: The quadratic factor in the small-frequency expansion of the conductivity 
function of the magnetic field for various values of the density. Left: / = 0. Right: / = 2. 



as a 



In figure 2, we show the behavior of quadratic term in relation to the magnetic field. 
We see that it behaves approximately as in the Simple Drude conductivity picture outlined 
in section B.l, with a few differences in the details. Essentially the second order terms 
for the diagonal and Hall conductivities start off at B — ► (In practice B = 10~ 5 ) at 
some negative value that is approximately proportional to the density and represents the 
relaxation time r -2 - where we notice the diverging relaxation time at po = = B that 
gave rise to the constant DC conductivity cr(u)) = j n (^'°) = £q = qiDTe or ct(uj) = De 
found in [28] and for a similar system first in [10]. At larger magnetic fields, it rises 
oc B 2 and becomes positive and then tails off after some maximum. The coefficient for the 
diagonal conductivity approaches a constant at large magnetic fields and the one of the 
Hall conductivity tails off approximately oc B~ l whilst the expectation from the Drude 
picture at constant r would have been oc B~ 2 - indicating at least a -B-dependence of 
relaxation time. The most striking feature is the "node" at which the correction term 
becomes independent of the density. In the Drude model, this would be the value of u c t at 
which the quadratic term vanishes. Looking in fig. 3 at how / shifts those curves, we find 




o"xx=o" w : p=0--- 
f=l,2,... — 
f=0.5,1.5,... — 



10. 12. B 14. 



O (7 

Figure 3: The quadratic factor in the small-frequency expansion of the conductivity as a 
function of the magnetic field for various values of /. Left: p — 0. Right: p = 2. 

that at small B, in the negative region in the case of finite p, that they are shifted towards 
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for increasing /, implying that the relaxation time increases, whereas for large values of 
B, they are shifted to larger values - which is simply an implication of the observation that 



UJ r T ~ 



Going a step further, we can check the generic predictions from section B.l. In fig. 5, 
we see that the ratio d| g 2 ^ ± at B ^ approaches precisely the prediction value | at large 
densities with a convergence rate that decreases with increasing / - even though we are in 
a completely different, i.e. strong coupling, regime. This also indicates that at large net 
densities and small /, we approach the classical Drude behavior, whereas for small densities 
or large /, we are in a completely different "phase". At large magnetic fields, however, this 
ratio does not become constant and depends significantly on /, but at least it seems that 
always ^Z- > . Furthermore, we can look at the location of the node, B cr it, which 



indicates the B value where 



«<7 



. .. = 0, i.e. where the peak turns into a minimum. For the 

B 2 B 2 

diagonal conductivity, we find that jfp « 0.342 at / = which converges to jfp « 0.397 
at large f. For the Hall conductivity, the value starts at 3.15, has a maximum of 3.90 

The variation in the ratio of those critical B 2 

B 



around / ~ 80 and then converges to 3.£ 

values is even smaller - between 9.18 and 9.84. If we were to associate tuj c 



as 



suggested from the DC conductivity in section 4.1.1, this is reasonably close to the values 
from the Drude model of |(u; c r) 2 and 3(lo c t) 2 . 




Figure 4: The ratio of the quadratic factors in the small-frequency expansion of the diagonal and 
Hall conductivities. Left: a ' CT " — °+ » - ln " 5 

9,V J 



at B = 10 5 as a function of p for various values of /. Right: 



. T n at large magnetic fields as a function of B for p = 10 and various values of /. The density 
p changes only the behavior at small magnetic fields and leaves the large- B tail unchanged. 



Looking in figure 4 at the quadratic term of cr xx at B 



0, where %^ becomes %^ 



—t , we find in fig. 4, that r is approximately proportional to the density, with a 
coefficient of r~ 2 f» 2.6/5. From another perspective, this means that the relaxation time 

— 1/2 

is approximately proportional to the mean distance between "quarks" , r oc p , but not 
the naive geometric mean free path in a system of weakly coupled particles. 

The proportionality coefficient is approximately r ~ 2.6ir^J^-. The /-dependence is 
not surprising, as increasing / appears to increase the relaxation time, which is consistent 
with a decreasing effective temperature that was a recurring theme in [28] . It is interesting 
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though that at large densities, the effect of / is only to shift the curves in fig. 4 and leaves 
the proportionality factor constant. 




Figure 6: The quadratic factor in the small-frequency expansion of the conductivity as a 
function of /. Top left: The coefficient of the diagonal conductivity for various values of the density 
and the magnetic field, right: For various values of the magnetic field at p — 1. Bottom left: The 
coefficient of the diagonal conductivity for various values of the density at B 6 {2,4}, right: The 
coefficient of the Hall conductivity for various values of the density at B S {0.5,4}. 



Looking at the coefficients as a func- 
tion of / in fig. 6 shows our observations 
from a different perspective. Essentially, 
the effect of / is to increase the relaxation 
time, and to decrease uj c t at fixed B. The 
most striking feature is the observation that 
we had above, that the coefficient in the 
Hall conductivity is proportional to / in 
regimes where it is negative, i.e. the "Hall 
peak" becomes narrower. 

Finally, we can look at the mass depen- 
dence in fig. 7. The biggest surprise from 
the Drude picture view is the quadratic de- 
pendence of second the expansion coeffi- 
cient on the mass. This indicates r oc rh~ 1 , which is somewhat counterintuitive since 
one would have thought that the relaxation time increases with increasing mass. If one 




Figure 5: The quadratic factor in the small- 
frequency expansion of the conductivity 



2<T 



as 



a function of the density for various values of /. 
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Figure 7: The quadratic factor in the small-frequency expansion of the conductivity as a 
function of to. Top left: The coefficient of the diagonal conductivity for various values of the 
density, right: For various values of /. Bottom left: The coefficient of the diagonal conductivity 
for various values of the B at / = 0.25, right: The coefficient of the Hall conductivity for various 
values of the magnetic field at p2. 

considers the Drude peak however to be a quasiparticle resonance, this is what one does 
classically expect since it means that the quasi particle becomes more stable at larger 
quark mass due to slower thermal motion and hence reduced collision rates. At vanishing 
density and different values of /, the result is also in contradiction with the free particle 
picture, since the DC conductivity is in a minimum at finite mass. There is an interesting 
maximum in the coefficient, which corresponds as / — > to the critical quark mass of the 
phase transition discussed in [31]. Hence, it occurs at the transition from the small-mass 
to the large-mass regime. This feature is even more apparent when plotting the coefficient 
against B for a small value of / = 0.25, where there is a small maximum around the critical 
mass. Looking at the Hall conductivity the regimes in B in which there is a Drude peak 
and in which there is a magnetoresistance minimum behave approximately like the pure 
Drude peak and magnetoresistance effects. It is an interesting curiosity, that the transition 
between those regimes receives a very small mass dependence. 

4.2 Large Temperatures: Diffusion limit 

In the diffusion limit, i.e. at uj <C q <C 1, we expect to be able to predict the transport prop- 
erties from the diffusion behavior, i.e. from the diffusion constant D and the susceptibility 
e because we expect the "mean free path" to be set by the temperature scale. 
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The diffusion constant was computed e.g. in [46] by studying the equations of motion 
of the gauge field in the gravity side to obtain Fick's Law, 



j(t, x) 



-DVj (t,x) 



(4.27) 



on the field theory side. The derivation in [46] is very instructive and can be followed also 
in our case in the presence of background fields. The expression for the diffusion constant 
is then slightly modified and yields 



D = — (-GV '-G tt G uu G xx ^j j 

7T 11 * 1 J q 



du 
^GG tt G uu 



(4.28) 



where we keep in mind that in our notation y/—G contains a factor of the u dependent 
coupling g~f f {u) = ^ f 2 + (1 - ^(u) 2 ) 2 . 

At M q = 0, this can be evaluated analytically and expressed in terms of hypergeometric 
functions as: 



(i + zVi + ^ + ^ + p 2 p 

ttTD = / du 

l + p + B 2 Jo 



1 + f + h 2 u ± 



(l + f)Jl + p + B 2 + ~p2 



2Fl 



l + P + B 2 
1 1 5 



1 + p + (5 2 + p 2 )n 4 ) \/l + {p + B 2 + p 2 )u* 

(4.29) 



4' 2' 4' 



(f z + B z + ~p z 



~P 2 



5.1 1.9. 

4' 2' ' 4' 



l;|;-(/ 2 + J B 2 + p 2 ), 



-B 2 



1+f 



5(1 + p) 



where 2F1 is the Gauss hypergeometric function, which is asymptotically in our case 
~ 4 ^^p + g2 + p2y/4 an d ^1 is an Appell hypergeometric function that is here at / = 

9IY5 /4} 2 

asymptotically ~ 5^(^2^2)5/4 . At F 7^ 0, the decay will be with a smaller, non-rational, 
power. In fig. 8, we see that the diffusion constant is for small / approximately propor- 




Figure 8: The diffusion constant ttT D as a function of the magnetic field and density for different 
values of / (left) and as a function of / and the density for different values of the magnetic field 
(right). 

tional to \fp whereas for large / the dependence is approximately linear. This may be due 
to the strong coupling because the usual classical geometric result for the diffusion constant 
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is proportional to the mean free path - which one expects to be inversely related to the 
density - and the mean free path should dominated by the baryon density at large baryon 
density, at least in weak-coupling intuition. However if we are for example in a superfluid, 
this intuition does obviously not apply anymore. 

At small /, the diffusion constant decays inversely proportional to the magnetic field, 
which represents the fact that charged particles in magnetic fields receive extra "drag" and 
become localized. At larger /, this decay slows down. Looking at the /-dependence, we see 
that the diffusion constant is approximately proportional to /, with an asymptotic slope 
that is independent from B. This contrasts to the dependence on p, which disappears at 
large B. 

Obtaining the permittivity is similarly straightforward. By definition [37, 38] 

e = lim C tt , (4.30) 

where it is understood that the limit oj — > is to be taken first. Taking the limit u>, q — > 
of the equation of motion for A t , (3.16), gives us a Poisson equation 



-GG G v 



l A' t (u))' = . (4.31) 



We note that this equation does not yield an appropriate infalling wave behavior near the 
horizon, but it is easy to see from the full equations for At and A x that for very small but 
finite wC^Cl, the behavior will be appropriately resolved near the horizon. Near the 
horizon, At and A x are strongly coupled, with At ~ ^A x , and A x follows an oscillatory 
behavior. To solve for A t , we then simply integrate (4.31) with A t = as a boundary 
condition at u — ► 1, which gives us readily the permittivity 



£ °A t 



-GG a G uu ^ 



e. n \ v = >u -° =: e e r , (4.32) 

fiduV=GG*G™ 



where (V—GG tt G uu ) u= o = 1. We can now see immediately, that the isotropic DC con- 
ductivity in section 4.1.1 is given by the diffusion result <j yy = a xx = eD and in the DC 
limit there is no contribution from other modes, as expected. However, in contrast to the 
remarkable result in the conformal case in [28, 10] where the conductivity was at all fre- 
quencies determined precisely by the diffusion behavior, the diffusion behavior is now only 
valid at small frequencies and receives corrections as we move away from oj = as outlined 
in section 4.1.2. Since the integral is the same as the one for the diffusion constant, we find 
that for M q = 0, we obtain the relative permittivity 



j? f 5.1 1.9. ( f2 , r>2 i ~2~i p z +B z 

pi 11 . 5 . + & + -(/ + B +p)i-^+r 

£ r ~ 2*1 T, , T, "(J +B + p J - p 



4> 2 '4' w ' ' r 7 5(1 + / 2 ) 

(4.33) 

In figure 9, we see the interesting fact that at large /, the relative permittivity becomes 
approximately constant. While one does not generically expect any specific dependence 
on the magnetic field, one would expect in a simple solid state model e r oc n and hence 
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Figure 9: The relative permittivity e r as a function of the magnetic field and density for different 
values of / (left) and as a function of / and the density for different values of the magnetic field 
(right). 



e r oc po at large p, which is realized here at large /, but at small / it is proportional to 
^fpo at large values of p. 

If we compare the results of this section to the Drude model reviewed in section B.l, 
we can identify from the Einstein relation (B.l) 

D 2 . 

(j, = — and n = irT eoe r ■ (4.34) 

4.3 T -> limit 

Next, let us look at the low temperature limit of q, Q S> 1. Here, we are interested in the 
equations near u = 0. The equations for A y and A x (3. 19), (3. 20) are identical in this limit, 
and become 

Ay + (Cj 2 -q 2 )A y = (4.35) 
as in the "conformal limit" in [28], up to order m 2 ^ or (e 2 , b 2 )^. The appropriate solution 



gives us the diagonal conductivity a yy = so\/l — q 2 /oj 2 and a xx 



4.3.1 Exponentially suppressed regime, q 2 S> 1 3> Co 2 , at small backgrounds, 

\Q\ > \P\,\B\ 

To study the low temperature limit more in detail, we start with the regime \p\, \B\ <C \q\ 
which is similar to the approximation in [28]. 

It is straightforward to see that the dominating term in the solution at finite h = 1 — u 2 
will still be A y ,A x ~ A° y ,A x e ±f i^- GXX / GU " . Using this to estimate the contribution of the 
cross-terms in (3. 19), (3. 20), we find that they are suppressed by a factor of q" 1 with respect 
to the dominant diagonal terms. In the near-horizon regime at u> 2 /q 2 </i<l, they are 
suppressed by a factor of h/q, and in the regime h <C io 2 /q 2 -C 1, they are suppressed by 
h 2 q/£j 2 . Hence, we can proceed as follows: First we will obtain the diagonal conductivity 
G yy {a xx follows similarly) by solving the homogeneous part, because the contribution from 
the cross-terms to the diagonal conductivity will be suppressed by the order of the square of 
the suppression of the cross-terms and can hence be safely ignored. Then we will compute 
the Hall conductivity from the inhomogeneous part. 
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Again, let us use the Ansatz A y = A®e$ which gives us 

c 2 + c + 



GGr"G«y)' {W=GG*(?»y) t G xx f G a 



1-QQuuQyy ' —Q QttQuuQxxQyy 

with the approximate result at h 3> co 2 /q 2 up to 0(1), £ 



co — q 



Quu \ Qxx 

(V^GG uu G yv )' 



= , (4.36) 



-GG uu Gyy 



Co := <7y gTHr, where we pick the negative sign corresponding to a solution that decays 

towards the horizon. Next, we take C 
remaining terms up to linear order in e 



(v^GG uu G yy )' c' 

towards the horizon. Next, we take ( = — Co — 2 V /= GG"" G^ 2<o 6 anc ^ S a * ner the 







6l'| + 2 C o 



^gQtuQxyy 



-G G a G uu G xx Gyy 



+ 



{^f^GG uu Gyy)' \ Co ( (V—GG uu G yv )' Co V 



+ 



v 2 V /3 GG UU G^ 2Co J 2Co 
=: e' — ea(u) — (5{u) . 



2J^GG uu Gyy 



+ 



2Co 



+ 



G u 



-to 



(4.37) 



The general solution to this equation is 

e = e /o" ^ eo + £ due" fo p(u) ) . 



(4.38) 



The second part is a small contribution £ M that is at most of order <j~ 2 , so we are only 
interested in the first part that evaluates to e = eoCoe 2 -^" d " Co( " ) , or e = e H Coe~ 2 ^ du ^ u \ 
en will be fixed in the region uj 2 /q 2 <C h <C 1, where there is an overlap between the 
asymptotic and near horizon solutions. 
At h <C 1, the equation becomes: 



-4^C + C 2 



4 A cD 2 



1 + + 

1 T 1j_f2 T" 



1+/ 2 T (l-*g) 2 +/ 2 



(4.39) 



As in [28] , this can be solved analytically in terms of hypergeometric functions and then 



be expanded for 



-fe > 1, giving us in the overlap region 



1+/2 (l_^2 )2+/ 2 



•Jh 



i+ 



1+/ 2 (l-*g) 2 + / 2 _|_ 



(4.40) 
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This solution connects nicely to the asymptotic region, even matching subleading terms in 
the overlap region, to give us 



Gvy 



-GG uu G yy )' 



Gyy 



G" 



-GG uu Gyy 



qvv 



+ 0(q 



G v 



Hence, the dissipative part of the diagonal conductivity reads to leading order 



Re a. 



yy 



(4.41) 



(4.42) 



where one could again interpret the result as having an "effective temperature" scale of 

T : 



1 , I G yy 

du\ 



T e ff vr J V G 
In the massless limit we can, as usual, find an analytic expression, which evaluates to 
T 

7V7 " VHX3/4) 



(4.43) 



4' T 4' 



(4.44) 



If we were to describe this result qualitatively as the behavior of a semiconductor, then 
the edge density of states would correspond to total density of baryons and anti-baryons in 
thermal equilibrium, and the difference -/V c — N v would correspond to the baryon density 




Figure 10: The "effective temperature" T e tt/T. Left: As a function of / and 
As a function of rh for various values of p and B. 



B 2 . Right: 



In fig. 10, we show how the effective temperature depends on the parameters of the 
defect. In addition to the dependence on /, that was previously found in [28], turning 
on a magnetic field or a finite density now raises the effective temperature approximately 
oc (-B 2 +p 2 ) 1 / 4 . Furthermore, we find that turning on a finite mass in some sense "enhances" 
the effect of the density and of / but the dependence on the mass in the presence of only 
B is not very significant. 
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Next, let us look at the off-diagonal terms. Do do so, we first need to write out the 
homogeneous part of the equation of motion for A x ='■ A x e-f *> x : 



, , (V^GG tt G xx )' 



G x 



-GG U G X 



which has the solution up to 0(1), ( x 



G^_ 

Quu \ Qxx 

(V=GG»G**)' & 



~2 ~2 
lo — q 



= , 



(4.45) 



2>/=SG**G"" ^ 2Co ± Where We again Pick the 

negative sign. The dominant terms in the full homogeneous solutions are then 



Ax 




and A h 



A 



y ^/V-GG uu Gyy( 



So Co 



(4.46) 



There are now two ways to determine the perturbative contribution coming from the cross 
terms. Either we again solve for the exponents ( - which would then contain factors of 
A® /A® - or we can directly take a pertubation for A y . Even though the latter one may 
seem most natural, in particular since the system is linear, we will use the first method 
since it gives us the result in a very neat way. Substituting (4.46) into the equation of 
motion for A y (3.19), we see that the equation for e (4.37) receives now an additional term 



(5{u) - (3(u) - Cx (V^GG tu G x y)' 



G uu A° x 
~G^A0 



(4.47) 



such that we obtain an extra contribution to e, taking only the leading term in (3 oc q 

-GG^ U G x ^ I 




G b 



G a 




(4.48) 



Now, if we look at the equations of motion (3. 19), (3. 20), we remind ourselves that in the 
near horizon geometry, the equations of motion for A y and A x look the same and the cross 
terms are suppressed by a factor of h/q with respect to the dominant terms. Hence the 
coupling occurs over the range u g]0, 1 — e] for small e and not in the near-horizon region. 
We keep e to regulate the asymptotic solution in the near horizon region, in which it is not 
valid. To capture the mixing then correctly, we fix e# at the horizon as above, such that 
we find 



Quu 
• thom. ~l" ^inh. 

where we absorbed a factor of 



o 2 ll £ Co 



G a G xx 



A y 



(4.49) 



G xx 



into €h- Because of the exponential factor in 



u=l— e 



the second term, the integral will be dominated around small u. The appropriate expansion 
gives us to leading order at « < 1 



£inh. 



A„ Ju 



qu 



U 



16g 4 (l + f 2 ) 2 



T(4,2qu) 



(4.50) 
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and in the limit u — ► 0, we find 

£5** = "I + | + "(A/ft + 0(1/,-) . (4.51) 

Finally, keeping A y \ u =o = A y fixed, we get the leading terms (ignoring the exponentially- 
suppressed terms) 

duAy = -^YTpqA° y - | P * A° x + 0{Ml/t) + 0(l/f) , (4.52) 

z (i + r) 7 <r 

and hence we can compute the Hall conductivity 
5A{ 
5A 



o yx = -e y/l + P — - 



= -eo I ~ - . 4 + 0{Mllf) + 0(1/ f) . (4.53) 



Following through the analysis attentively, one can also see that the imaginary part of the 
Hall conductivity is exponentially suppressed by a factor of e~ q / Te ff . 

This result is remarkable, since the diagonal part of the dissipative conductivity is 
heavily suppressed with a factor e~ q / Te -ff , while the off-diagonal part is only suppressed by 
a factor of T 4 /g 4 . This reflects the fact that at small temperatures we approach conformal 
symmetry in the field theory and hence the form of the conductivity in [44] that was 
discussed in [28]. Having a purely off-diagonal conductivity is not surprising as it is for 
example the case on the Hall plateaus in the quantum Hall effect or as we demonstrated 
above in intrinsic semiconductors at small temperatures. 

Interestingly it occurs also in intrinsic semiconductors at low temperatures. There the 
absence of defects and the highly suppressed charge carrier density cause the relaxation 
time r to diverge, whilst the carrier mobility remains approximately unchanged. Hence 
the diagonal conductivity (B.2) is suppressed, while the factor uj c t in the Hall conductivity 
(B.3) causes the Hall conductivity to remain finite. 

4.3.2 Exponentially suppressed regime at large backgrounds, q 2 » |/5|,|-B| » 
\q\ > 1 

This regime is slightly more non-trivial, because now £o splits into three regimes (for 
simplicity at M q = 0): 

u<-^=,^: Co ~ 9\A + P , (4.54) 



B 



1 ! ^<n«l: Co - , q ^ l + lL^- and (4.55) 



h < 1 : Co ~ QV 1 + P ^_ _ (456) 



y/p+B 2 +p*h 



In the asymptotic region, the solution is dominated by the decaying exponential, whereas 
in the near horizon region, it can be written in terms of the coordinate s from section 4.1.1 
as 

A y = A° y e ius = A°e uaS ° v=gg™gvv , (4.57) 
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, _ /rT - ? 2 > /g 2 (/ 2 + (l-^) 2 ) + (l+/ 2 )(p 2 +/ ii + (l-^) 2 ) . . .. , 

where v = wyl + j z l+ ^2 + p • t ms solution were to overlap 

with the "tail" of the asymptotic solution, u 4 > f 2 + p 2 + B 2 , we could match them at 
some 1 » «f > (/ 2 + p 2 + B 2 ) 1 / 4 . The fact that they do not overlap, however, can be 



simply seen from the different u dependence of d u s = ^—^^ uu( ^ an d Co = Qy gSr- Hence, 
whatever we try now, the conductivity will disagree by some finite factor. One way to 
pretend that they do overlap is to simply set h = 1 — u 4 — > 1 and u — > 1 , which corresponds 
to extending the -4= , 4= < u <C 1 region towards the horizon and the near-horizon limit 

into the intermediate region. Matching the solutions under these conditions and ignoring 
the second part of the solution for e (4.38) gives us 



1 gvv 



* = -&\ \ V GUU (4-58) 



1 + f 2 V G 

with the corresponding conductivity 



1 + f 2 + p 2 + B 2 _ 2f}J ;; ./• 



^ = -* £ o| + • . . + e y 1+ ^ 2 e V g uu + . . . . (4.59) 

Certainly we cannot trust any 0(1) and polynomial factors, but the point to make here 
is that we should still expect the exponential suppression from the effective temperature 
(4.43). This is because the integral in the exponent is dominated by the region in which 
Co does indeed dominate the solution, we took care of the deep near horizon region and 
elsewhere there are no "large" terms in the equations of motion. 

In principle one can also try to find a solution in the regime ^ ~ u ~ 1 an< ^ 
then "glue" it to the near horizon and asymptotic solutions to gain more accuracy, but 
there is limited insight to be learned from this and it would be much more tedious than 
the calculations in 4.3.1. 

The Hall conductivity will still be dominated by the asymptotic regime since the mixing 
from the near horizon region is exponentially suppressed. In principle, it is then still 
&yx = — £0 1 ( 1+ jf^s/2g4 1 however we need to note that the integral that was computed in 

section 4.3.1 is dominated around the maximum of u 3 e ^ 2 ^V 1+ f 2u at u max = j and 

decays then also on the scale 5u = 2 -^/^j^ • C an d hence also the exponential suppression 

however start to change around u s t op ~ (/ 2 +/> 2 +i? 2 )~ 1 / 4 . To properly evaluate the integral 
in the approximation (4.50), u s top need to be significantly larger than u max - otherwise the 
"tail" of the polynomial term in the integral will not be sufficiently suppressed and will give 
a finite contribution to the result, which will greatly overestimate the Hall conductivity. 
Certainly one can always use the full integral 



(4.60) 



but this is certainly a somewhat less insightful result and cannot be computed analytically. 
In figure 11, we demonstrate the boundaries between the different regimes. To see how 
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Figure 11: Left: The check to the deviation from the T — > limit as described in the text. Right: 
The numerical estimate j^jj <~ — 7r<9,jln ( compared to the analytical result. 



far the regime of the previous section reaches, we plot ln(a yy ) — In (cr y y Sym '^ against 

fiff ■ <Jyy Sym "^ is defined taking In — ~\ to be the linear expansion of In (-f 1 -^ at 

large values of q approximately where a yy ~ 10~ 15 £0) shortly before the numerics fail. 

= n corresponds to q = p n . This is sensitive to changes in the factor in front of the 
suppressed conductivities at large q. To check for the overall limit of the exponentially 
suppressed regime, we look directly at a yy — a y a y yrri "\ It is easy to see that the boundaries 
at approximately p, B ~ q and p,B^q 2 are verified. Using that data, we also looked at 
the Oil) factor ir in front of the exponential term in the conductivity in (4.42). It turns 
out that for our values of tj^— ~ 35, the numerically estimated factor varied from ~ 3.8 
at p = to 15 at p = 140, where the numerics carried us only up to q ~ p. At the latter 
values, we could not expect close agreement because we were outside the regime that we 
considered in section 4.3.1 - and the value of 3.8 seems reasonably close to ir. 

4.3.3 Dominantly large backgrounds p, B » q 2 S> 1 

As p, B > q 2 , we find that J* ( ~ 2g(/ 2 + p 2 + B 2 )^ 1 ^ <C 1, and hence the exponential 
suppression factor disappears. Furthermore, as indicated above, our estimate for the Hall 
conductivity does not apply anymore, as it will be dominated by the region of u in which 
the assumption (o 3> 1 does not apply anymore. 

Looking at the problem in another way in terms of the coordinate s from section 4.1.1 
for the DC conductivity with the equation for A y (4.1) gives us the relevant term 

d u d s A y = ... + V^GG yy G xx q 2 A y . (4.61) 



Integrating this analytically in the massless case and for for sufficiently slowly varying 

(i+/ 2 )g 2 ' 

(/ 2 +B 2 +p 2 )l/4 



d s A y < Ay, i.e. v < 1, we find 6 d s A y J* du d s A y < TjStir^Tjl G ^. A second 



integration gives < — ^=gH= and for larger frequencies u, the oscillatory behavior 
of A y implies that the integral is further suppressed by a factor 1/uj. 
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Hence, for p,B^> q 2 , the real part of the conductivity is dominated by the results 
of the isotropic case in section 4.1. The inductive (imaginary) part of the conductivity 
obviously still receives the term a y y ~ — i£o^- 

5. Numerical Results 

In this section, we study information that can be derived from computing the correlators 
numerically, in particular the overview of the frequency-dependent conductivity (i.e. the 
spectral curves), the diffusion and relaxation behavior in the hydrodynamic regime, and 
the spectrum of quasi-particles. 

5.1 Spectral Curves 

In this section, we present the conductivity spectrum in the presence of various background 
quantities. 




Figure 12: The real, diagonal part of the isotropic conductivity at q = as a function of frequency 
and net baryon density (top) or magnetic field (bottom) for / — (left) and / — 2 (right). 

First, let us look at the case of finite density and magnetic field alone in fig. 12. From 
the result of the electromagnetic duality in section 3.1, and the very generic results for the 
hall conductivity in section B.3, we expect to see a sequence of resonances, in which maxima 
and minima are exchanged between the case of finite net baryon density and magnetic field. 
The fact that "plasmon" (finite-density) resonances are relatively strong is not surprising 



- 33 - 



since this is a strongly coupled system - and plasmons are a finite coupling effect. The 
small-frequency regime reflects very well the classical Drude model expectations and the 
small- frequency expansion from section 4.1 - with the Drude peak and magnetoresistance. 
Looking at the resonances, we find that they are approximately equally spaced at nco p or 
(n — l/2)u} c , respectively - and they decay quickly. Comparing this to what we learned in 
section B.3 reveals interesting information about the quasi-particles that carry the current: 
They a) must be massive and b) do not consist of chiral fermions, in sharp contrast to 
graphene [41]. It is also interesting to see that the amplitude seems to be again decaying 
exponentially as in the resonance on the width of the defect that were studied in [28] . The 
frequencies uj c and u p are, however, not proportional to B and, respectively, p and even 
start off at a finite value. In terms of generic weak-coupling intuition, this would need to be 
explained by a non-linear magnetization behavior and non-linear chemical potential, and 
in terms of a changing mass of the quasi-particles that carry the current. 




Figure 13: The real, diagonal part of the isotropic conductivity at q — as a function of frequency 
and magnetic field for B = 4 (left) and p = 4 (right). 

Looking at the /-dependence in figure 13, we find that on the one hand, increasing /, 
i.e. an increasing width z max or stronger "confining potential", reduces the amplitude of 
the resonances at small frequencies. This is consistent with the value that we found for the 
DC conductivity (4.1.1). In contrast to this, we find that the suppression of the resonances 
with increasing frequencies decreases with increasing / and we can see the tower of modes, 
that is at small amplitudes hinted at by the a = e$ lines in the plot. This agrees with 
the effective temperature (4.43) that decreases proportionally to / _1 ^ 2 . Secondly, we find 
that the parameters lo c and u p decrease with increasing /, which we can again explain by 
a non-linear behavior of the response functions or by an /-dependent quasiparticle mass. 

To see what happens when we turn on a finite wavenumber of the perturbations, we 
look at fig. 14, where we show a few of the higher resonances at / = 4 in order to see how 
they depend on the wavenumber q. Looking at the plots on the left, there seems to be only 
a small difference between the behavior of the Landau levels and plasmons. This difference 
becomes however very significant when one plots the "normalized" conductivity, , a 
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Figure 14: The real part of the conductivity a yy at varying q, f = 2 and p = 4 (top) or B = 4 

j (left) and rescalec 

frequency \JCj 2 — q 2 (right). 



(bottom) as a function of u> (left) and rescaled as , <Tyy as a function of the "rest-frame" 



as a function of the "rest-frame frequency" -\J 6j 2 — q 2 as it was done in [28]. Then, we see 
that the density resonances connect smoothly to the resonances in the optical regime (i.e. 
above the conduction threshold u = q) in the "semiconductor" case at q S> 1. Certainly the 
statement about the continuity of the pole or "resonance" at -^/tD 2 — q 2 = is somewhat 
meaningless, since this arises always due to the rescaling (at finite temperatures), but only 
says that the correlator is finite at u = q and does not imply a pole in the correlator. 
The magnetic resonances, however, seem to be discontinuous - the n = Landau level 
seems to disappear, when the \/oj 2 — q 2 = pole arises, and the higher resonances behave 
in a non-monotonic way. In order to see more in detail where this discontinuity comes 
from, we can look at the B and /j-dependence at a finite wavenumber q = ir/2 and finite 
/ = 4 in figure 15. There we see that we start off with the "bare defect" and its finite- 
width resonances, and as we turn on the net baryon density, they shift smoothly, as if we 
were to decrease the width of the defect. As we turn on a magnetic field, while there is 
still no apparent splitting of resonances - as one might expect if new kinds of resonances 
are turned on, they are not monotonically connected. This implies that there are some 
non-monotonous changes in the residue and location of the poles. It can be easily seen 
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Figure 15: 



Re a v 



at q = ir/2 and / = 4 as a function of the "rest-frame" frequency \Ju)' 2 — q 2 



. Left: As a function of the magnetic field. Right: As a function of the density. 



from the electromagnetic duality in the plain defect, that at wavenumber q = 0, there can 
be only one pole, which is at u) = 0, and hence, assuming continuity, the residue of the 
poles from the finite-width resonances must be proportional to q. On the gravity side, this 
corresponds to the fact that at q = there is only one mode function in the gauge field 
and the equations for A y and A x arethesame, but at finite q, the equations for A y (3.19) 
and A x (3.20) become different. Hence we find two distinct mode functions. The same 
argument applies for turning on p or B. This also reflects the fact that generically, the 
density of states of Landau levels (B.14) is proportional to the magnetic field. 

A rough, argument in the field theory is that turning on q corresponds to introducing 
an inhomogeneity in the x direction. Hence, the U{\) perturbations become localized 
in that direction, whereas they are not localized in the y direction. Plasmons are not 
generically localized, so they do not change this configuration. Landau levels however are 
intrinsically localized quasi-particles, so they break translation invariance also in the y 
direction and change the pattern of the resonances less smoothly. From EM duality, we 
know that magnetic resonances in a xx connect smoothly to the finite- width resonances and 
the density resonances connect less smoothly. This is precisely because in the x direction 
the translation invariance of the plasmons becomes broken by finite q, whereas the Landau 
levels were already localized. 

After studying the effects of having either B or p turned on, let us look at the case 
when they appear simultaneously in fig. 16. In the diagonal part of the conductivity tensor, 
we see that the magnetic- or density resonances split in two as we turn on a net density 
or magnetic field, respectively. It is interesting, that there is no "tower" of excitations 
splitting off each resonance and that the mean frequency of each "split level" charges 
only by a small amount. Furthermore, we find that at each resonance in the diagonal 
conductivity, the Hall conductivity changes sign, at least for the first two resonances. This 
is just the continuation of what one expects classically for the first resonance as we saw 
in section B.l. It is also what one expects semi-classically, if the split states have either 
positive or negative magnetic moment, carrying a total net Hall current similar to the edge 
current in the Quantum Hall effect. By continuity, this implies that the plasmons (at zero 
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Figure 16: The real part of the isotropic conductivity at q = 0, B = 4 and / = 2 as a function 
of frequency and density (top) and at p = 4 as a function of the magnetic field (bottom). Left: 
diagonal part of the conductivity tensor. Right: Hall conductivity. 

magnetic field) and the Landau levels (at zero net "charge" density) have vanishing net 
magnetic moment and equal degeneracy (2). It is also worth noting that it is impossible 
to have any of the resonances cross uj = no matter how much one tunes the parameters, 
which clearly indicates that the system has no Fermi level in the classical sense. Otherwise 
we would see Landau levels crossing u = 0. 

For completeness, we can look at the /-dependence of the Hall effect in figure 17. This 
confirms our observations of the relation between the resonances in the diagonal part of 
the conductivity and the Hall conductivity. In the regime of highly suppressed resonances 
this appears through their periodicity that differ by a factor of 2. We can also see that 
their frequencies roughly behave as the resonance frequencies of the plasmons and Landau 
levels. 

Finally, let us look at the mass dependence. In fig. 18, we look at the conductivity at 
q = p = B = 0, where we actually see the DC conductivity from (4.1.1). At / = 0, we see 
a significant change of the conductivity with a resonance around w ~ 1.7 as m approaches 
the critical mass of the phase transition. This m-dependence is suppressed at finite /, and 
at / = 2, the most significant change takes place only over rh ~ ... 4 - simply because 
it depends roughly on f 2 + (1 — \I/ 2 ) 2 and not on the mass directly, such that the mass 
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Figure 17: The real part of the isotropic conductivity at q — 0, B — 4 and p — 4 as a function of 
frequency and /. Left: diagonal part of the conductivity tensor. Right: Hall conductivity. 




Figure 18: The real, diagonal part of the isotropic conductivity at q — as a function of frequency 
and quark mass for / = 0, rh G [0, 1.18] (left) and / = 0, m G [0, 8] (right). 

dependence becomes "frozen" as \&o becomes close to 1. In contrast to this indirect mass 
dependence, the location of the very shallow maximum seems to be roughly proportional to 
rh. This gives some nice insight into the IR and UV dependence of the underlying physics. 
Processes that take place at small energies, i.e. in the IR will be dominated by gravity 
background near u = 1 , and hence depend on \&o an d show most of their mass dependence 
in the regime of rh ~ 0{1). Effects that depend on high energies, i.e. the UV, however 
depend on the background near u = and hence depend on rh (and only to subleading 
order) on c. We can observe the influence of the quark mass on the fmite-t/ resonances in 
figure 19. There we see that the gap between the resonances is roughly proportional to 
rh~ x at large rh and the change starts oc m 2 at small rh - as one does generically expect for 
a relativistic system. As naively expected, the resonances are also narrower at large rh and 
their amplitude increases. If we look at the overall level of the conductivity (i.e. ignore the 
resonances) there seems to be the correction that we found at q = 0, now as a correction 
to the background around which the resonances take place at small yju 2 — q 2 . This also 
agrees with the picture that we see at / = 0. Looking at the un-scaled conductivity, 
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Figure 19: Top: The real part of the normalized 
conductivity , <Jyy = at q = ir/2 as a function 

of "rest-frame" frequency \/uj 2 — q 2 and quark 
mass at q = vr/2 for / = 0, m G [0, 1.18] (left) 
and / = 2, to G [0,8] (right). Bottom: CT yy for 
/ = 2, tog [0,8]. 



at the bottom in fig. 19, we also see that the conductivity approaches the t — » limit, 
= £oRe a/1 — q 2 /u> 2 as we increase the quark mass. 




Figure 20: The real, diagonal part of the isotropic conductivity at q = as a function of frequency 
and quark mass at / = 2 for B = 4 (left) and p — 4 (right). 

These generic effects of turning on M q can also be seen in the plasmons and Landau 
levels, and in the Hall effect, in figures 20 and 21. Again, we see that on the one hand, the 
resonances become more stable at large m, and on the other hand that the energy levels 
receive at small m a small correction oc rh 2 and at large masses scale oc m _1 , just like u> c 
and uj p do classically. 
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Figure 21: The real, diagonal (left) and Hall (right) part of the isotropic conductivity at q = 
and / = 2. Top: As a function of rh at B — 4, p — 4. Bottom: As a function of p at -B = 4 and 
m = 8. 
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Figure 22: a yy — a yy and a xy at q = 0, rh = 32, 
p = 32 and B = 24. 



Finally, we can turn on a large mass 
(in this case rh = 32) in order to study the 
structure of the Hall effect more rigorously. 
In fig. 22, we see that the Hall conductiv- 
ity has a small overall positive (or negative 
if we rather look at a yx or negative Bp) 
background, and there are poles with alter- 
nating residue, each precisely located at a 
maximum of the diagonal part of the con- 
ductivity. This supports exactly our sug- 
gestion above that the Hall current is car- 
ried collectively by localized states with net 
positive or negative magnetic moment. 
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5.2 Small frequency regime 

In this section, we look at the behavior of the purely dissipative poles of the correla- 
tor C xx on the imaginary axis, that dominate the conductivity at small frequencies and 
wavenumbers. Our particular interest is how they influence the DC conductivity and how 
the transition to "semiconductor-like" in the quasiparticle regime at larger wavenumbers 
occurs, i.e. how those poles disappear. 

The numerical strategy behind locating the poles is reasonably straightforward. First, 
we divide the imaginary axis in three regions, based on an educated guess, and localize the 
poles in these regions in a simple recursive process at some initial wavenumber, magnetic 
field, quark mass and density. Then, we can identify regions around those poles that 
allow us to "track" them as we change the parameters, without having to scan the whole 
imaginary axis. One caveat though is, that it is numerically increasingly difficult to find the 
poles as their residue decreases, so we keep a minimum wavenumber (we will use q > 0.02) 
to always find the "middle" pole. We may also "lose track" of poles if their residue becomes 
too small. The other caveat is that with our rudimentary method, we need to filter the 
result afterwards for whether a suspected pole is a pole or just a local extremum or noise. 
In most cases the distinction is obvious, but in some cases we will look at the value of 
the residue that we estimate. Furthermore, since this process is reasonably numerically 
intensive, we will limit the computations to a few examples. 

Applying accurate numerics, we find that there are only at most three purely dissipative 
poles; the diffusion pole at small imaginary values of a) <C 1, and two more rapidly decaying 
poles at cDO(l). We recall that in [28], there were found the first two of those three poles 
- the diffusion pole and one corresponding to decay on thermal scales. Those poles were 
found to move along the imaginary axis as we increase the wavenumber, meet at some 
critical wavenumber, and for more short scale excitations turn into massive quasi-particles. 
Obviously, at q = B = p = 0, there can only be the diffusion pole because then the 
electromagnetic duality together with isotropy restricts £q 1 a £ { — 1, 1}. As we tune these 
quantities to zero, the residue of both poles vanishes. One of them just disappears to a 
constant conductivity, while the other one turns into a unit step function of frequency in 
the conductivity. 

First, let us look at the poles in the presence of a finite density. In fig. 23, we show the 
frequency of the diffusion pole and the second pole as a function of p, and we see that the 
lower pole follows, at small q, the diffusion behavior, and then, at some critical wavenumber, 
they merge and we have again the branch cut with the transition from the dissipative to 
quasiparticle behavior. It is interesting to note that even beyond the diffusion behavior, the 
curves agree reasonably closely upon the appropriate rescaling with the diffusion constant. 
In fig. 23, we also plot the critical wavenumber as a function of p. We see that there is 
no length scale in our system that fits it particularly well compared to a simple 7r _1 /5 -1 / 2 
approximation - even though -i seems to fit best asymptotically. As in the case of the 
diffusion constant, we note that the critical wavenumber is proportional to l/y/~p, and ^fp 
is approximately the mean separation of the quarks. So at a small net quark density, q c 
is dominated by scattering off gluons and quarks from the thermal equilibrium, and at a 
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Figure 23: Left: Imaginary frequency of the diffusion pole and the first higher pole as a function 
of q for different values of p. The frequency and wavenumber are scaled with D, such that the 
diffusion equation is Dlu = —i(Dq) 2 . The dot indicates the point where there is the branch cut, 
and the mode frequency of the poles gain a real part to become quasi-particles. Right: The critical 
wavenumber as a function of p, compared to various length scales of the problem: the effective 
temperature, diffusion constant and the electric permittivity. 
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large net quark density, it is the baryon density that sets this scale. 

The imaginary frequency of the third pole is slightly increasing with increasing q, but 
it has a very small residue that decreases with increasing q. Hence it can only be seen at 
q ~ O(0. 1) and we do not plot it here. 

To study the nature of the second pole, 
we look in figure 24 at how the location of 
that pole depends on the net baryon den- 
sity. We find that this too is proportional 
to p~ 1 / 2 at large p, and approaches some fi- 
nite value at small p. Also, the dependence 
on / dominates only the small- p regime. In 
that figure, we also compare the location 
of this pole to the inverse of the relaxation 
time r _1 that we obtain from 




Figure 24: The location of the second pole on 
the imaginary axis as a function of the net baryon 
density for various values of /. For comparison, 



as defined in section B.l and computed in 
4.1.2. We find that for large p, they are in 
perfect agreement, whereas for small p, the 
inverse relaxation time diverges. For clar- 
ity, we show only the relaxation time for 
/ = 0. The relaxation time for larger / is 
in worse agreement (r _1 becomes larger at 
small p), however the agreement at large p is equally good. This disagreement at small p 
reflects the special conformal nature of the system at p = 0, with the constant conductivity 
from the electromagnetic duality in [28]. This causes the relaxation time that we com- 
puted from the conductivity to vanish, while we can obviously expect that any excitation 
still decays on a finite timescale as dictated for example by causality. The reason why 



we show the inverse correlation time t _1 com- 
puted from the conductivity for / = only. The 
other curves for the relaxation time will be in 
worse agreement at small frequencies. 
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we do not see this relaxation time in the conductivity at q = p = and hence why the 
constant conductivity does not violate causality is that the residue of this relaxation pole 
vanishes in the isotropic limit at vanishing density Furthermore, we notice that this theme 
of exact convergence to the Drude model ((B.2) has precisely a pole at w = —it' 1 ) at 
large p is recurrent and was already seen in figure 5 in section 4.1.2. From naive intuition 
about weak coupling, one might be puzzled as to why the relaxation time scales as y^p, 
i.e. proportional to the inverse of mean separation between the quarks whereas in a simple 
geometric weakly coupled model, the relaxation time is proportional to the mean free path, 
that is proportional to the density. Because this system is strongly coupled and there are 
long-range correlations, however, this intuition breaks down. 




Figure 25: Left: The location of the second relaxation pole as a function of the net baryon density 
for various values of /. Right: The logarithm of the residue of the second relaxation pole. 

This case of very small q = 0.02, varying p and / is also a good example to study 
what happens to the third pole. To do that, we can look at fig. 25, which shows the 
location of this higher relaxation pole and (the logarithm of) its residue as a function of p, 
at various values of /. We see that the residue decreases exponentially with increasing p, 
while the poles move to higher imaginary frequencies, until they can't be tracked anymore. 
To illustrate this better, the frequency in the region of p in which there is no reliable residue 
information anymore is plotted dashed. This is because the accuracy of the location of the 
pole that is necessary to determine the residue with our methods is of the order of the 
residue. Obviously, there is no guarantee that the extremum that one finds actually is a 
pole if there is not good enough data to compute the residue. One could argue that we 
could go to much higher accuracy since the accuracy grows exponentially with the number 
of steps. However the number of steps needed to track the poles grows with the inverse of 
the residue, i.e. exponentially with increasing p or q. If we look at the density dependence 
of this pole, we find that its imaginary frequency increases with increasing density - which is 
just what we expect for a naive model of weakly-coupled particles with finite cross sections 
a. In particular, for small cross sections in d dimensions, cr 1 /^- 1 ) <C p^ 1 ^, one expects 
classically r^ 1 ~ vap, where v is some characteristic speed. Hence for larger densities or 
larger cross sections, the scaling would naively approach yj~p. This is just what we see in 
fig. 25. Hence, despite the limited accuracy and reliability in tracking these poles, we can 
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safely associate this pole with a classical, weakly-coupled relaxation behavior. 




Figure 26: Left: The location of the relaxation poles as a function of the magnetic field for various 
values of /. Right: For different wavenumbers. 

A behavior very much in contrast to the case of turning on the finite density can be 
found in in the left of fig. 26, where we show how the relaxation poles behave at small 
wavenumber in the presence of a magnetic field. We see that with increasing magnetic field, 
the relaxation poles merge at some critical magnetic field, B c , and then turn into the first 
Landau level. This contrasts to the classical Drude-model analysis, where the magnetic pole 
moves away from the imaginary axis as soon as the magnetic field is turned on. Essentially 
what is happening is that the creation of the first Landau level is inhibited below B c 
because of the strong coupling. If we assume a crude model, in which the frequencies are 

given by -i cL 2 ± V^c - — 4 . then we see the reason for the dependence of 
the curves on f. Now, remember that we found in section 4.1.1, that uj c t ~ , B , such 

that B c is approximately given by the ratio of the relaxation times at vanishing magnetic 
x /T+f2 f r- 1 \ 

field, B c ~ -*-2 — I ;pr — 1 ) • This ratio depends non-monotonically on /, because 

apparently the location of the classical relaxation pole is not closely related to the location 
of the first relaxation pole. Obviously, this generic behavior is not exact, but provides a 
rather qualitative description. 

This magnetic dependence is consistent with the dependence of the relaxation poles 
on the wavenumber that we show in the right of fig. 26. In fact, now it is most apparent 
that the behavior is reasonably well-described in terms of t c \Jt only. This can be seen 

by computing -i- ( ^ — 1 J . For the data in fig. 26, this ratio is approximately 
c V T J B=0 

constant, ranging from 4.2 at q = 0.02 to 5.3 at q = 1.2. Certainly it is not close to 1, but 
we could not expect this, as we chose our expression only as an example of how a branch 
cut in the solution for the location of the poles as a function of the background parameters 
can look like. This shows us however how suppressing the first Landau level is related to 
strong coupling, as the critical magnetic field is proportional in some approximation to the 
difference between the classical and strong-coupling inverse relaxation times. 

Finally, we can look in fig. 27 at the dependence of the location of all three poles 
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Figure 27: The location of the poles on the imaginary axis as a function of the wavenumber for 
different values of the magnetic fields. The frequency and wavenumber are scaled with the diffusion 
constant, as the diffusion behavior can be written as iDuj = (Dq) 2 . 

on the wavenumber at varying magnetic field. At vanishing magnetic field, we start off 
with the system in which there is the hydrodynamic to quasiparticle transition, and there 
exists always the classical relaxation pole, that moves towards larger r^ 1 with increasing 
q. Going beyond the critical magnetic field, we see that the relaxation poles re-appear at 
some wavenumber qs < q c that increases with increasing magnetic field. Beyond some 
second critical magnetic field, at which qs = q c , there is only one imaginary pole, that 
starts off at small q as a diffusion pole and turns at large q into the classical relaxation 
pole. Interestingly, at large wavenumbers, we always see only the classical relaxation pole, 
so the "strong" relaxation pole is an effect that arises only at small wavenumbers, i.e. 
long-distance perturbations, which is what one actually expects, because of the diverging 
correlation length. The effect that the magnetic field inhibits relaxation on large length 
scales is precisely what we expect because of the localizing effect of the magnetic field 
on charged particles. This behavior between the critical magnetic fields, in which a pair 
quasiparticle poles (with positive and negative real part of the frequency) turns into a pair 
of relaxation poles and then into a different quasiparticle pole nicely reflects the transition 
between the regime dominated by Landau levels and the regime dominated by resonances on 
the width of the defect, that we observed in the previous section. Furthermore, we see that 
under the appropriate rescaling with the diffusion constant, the curves agree even beyond 
the actual behavior proportional to q 2 - indicating that q c is reasonably well-described by 
the lengthscale from the diffusion constant. The attentive reader will notice the "hole" in 
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the plot near q c . This arises because it is difficult to track the poles in this regime using 
our method of scanning the frequency for each value of q in this case. Obviously, swapping 
that order should allow us to determine the location of the poles also in this regime. 
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Figure 28: The imaginary frequencies of the relaxation poles as a function of the quark mass. 
Left: For various values of /. Right: For various values of q. 



The effect of turning on a finite quark mass is very similar to the case of the background 
magnetic field, as we see in fig. 28. The most significant difference is in the /-dependence, 
as the critical mass at which the relaxation frequencies receive a real part increases quickly 
as we increase /. This is due to the fact that the hydrodynamics is dominated by physics 
in the IR, i.e. at small radii and hence depends on rather than the mass directly and 
the quark mass as a function of 'I'o increases with increasing /. 

Looking at the picture of the location 
of the poles as a function of the wavenum- 
ber in fig. 29, we find again that the system 
at small imaginary frequencies and small 
wavenumbers is well-described by the dif- 
fusion behavior, however not up to as high 
wave numbers as in the case of the back- 
ground magnetic field. The critical mass 
above which there is always only the classi- 
cal relaxation pole is just below rh ~ 0.75. 

5.3 Landau Levels and Plasmons 




Dqo.8 



Figure 29: The location of the poles on the imag- 
inary axis as a function of the wavenumber for 
different values of the quark mass. The frequency 
and wavenumber are scaled with the diffusion con- 
stant. For comparison, the critical mass for the 
embedding is fa <~ 1.196. 



In this subsection, we will try to shed some 
light on the nature of the magnetic and 
density resonances that we observed in sec- 
tion 5.1, by studying the corresponding quasi- 
particle poles. For simplicity, we will focus 
on the isotropic (q = 0) case. In principle, there are again different methods of estimating 
the location of the poles. The least reliable method is simply fitting Lorentzians to the 
resonances. However it cannot give the right answer, as we expect a sequence of infinitely 
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many poles with separation uq — 270, and if we consider the n th resonance and provided vq 
and 70 are of the same order, we need to consider more than 0(n) neighboring poles. The 
more precise method involving only the data on the real (frequency) axis is then as in [28] 
based on assuming an appropriate sequence of poles, summing it (ideally analytically), and 
fitting the parameters locally around each maximum - assuming they vary slowly enough, 
such that the "backreaction" from the varying parameters is sufficiently suppressed. The 
third method is simply trying to fit the poles by scanning an area in the complex frequency 
plane using an appropriate guess obtained from the data on the real axis. Then, we can 
use the usual recursion to find the poles. The most time consuming step in that method 
is to scan the search area for the first time, since we can not assume that the poles are the 
only local extrema. For the former technique, we use the Ansatz 

^1/ n(P + «7o) 2 ,~ . ._ s . n(£ -«7o) 2 ... .-A n 

that was also used in [28]. This Ansatz basically says that all the poles are located with 
equal spacing on a straight line at v n = n(±vo — «7o) (or v n = (n — 1/2) (±^o — *7o) for the 
Landau levels), with residues (or f 1 n "i/ 2 ) • 

However, it turns out that methods to try to fit the resulting analytic expression 

/ ~ x sinlg% - isinhlg% 

Im C vv = e Im Pcot vr = e Im v > ( 5 - 2 ) 

VV V vo + i%) cosh - cos & V ' 

sinhlg% 

Re or™ = e 9 - — 5 - , (5.3) 

ra cD C o S h IS - cos & V ' 

^0+7o ^o+7o 

are unreliable in this case, mainly because of the periodicity of the expression. 

From (5.3), we can then obtain the parameters for the poles from the n th local max- 
imum, a n , Lo n and its neighboring minima cJ n _ 1 / 2) a n+\/2- Assuming a n _i/ 2 = & n +i/2 = 
a m in, the exact result is in the case of plasmons 

--o £o u n = JtpL (5 .4) 



sinh 2 "-l f o 

"o+7o 

and similar for Landau levels. This expression can be trivially inverted to obtain vq 
and 7o- Taking into account that cr n _ l / 2 / o" n +i/2) the correction when using o" m j ra = 

2-2/ ~\ — 1 

\ (< T n-i/2 + a n+i/2) wm be of order ( cosh ^ ) , i.e. it will only be significant for 



'n-l/2 T" <J n +l/2 

the first few poles. However this seems to be a less bothersome shortcoming to take than 
"misfitting" poles because of the periodicity. Note that the two terms in the error term, 
make sure that beyond n = 1 always either of them gives us a good suppression. 

In fig. 30, we compare the two methods to obtain the poles - the "exact" direct 
search and the "approximate" result from fitting to the Ansatz - for various values of /, 
B and p. In the direct search, we used the maximum value of the spectral function in the 
last recursion as an indicator as to whether the pole has been found, with a threshold at 
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Figure 30: The location of the poles u> n — v n — ij n in the complex frequency plane from the 
different estimates as described in the text. For orientation purposes, the lines connect the poles of 
different order in the same background. Left: Finite p = 2, Right: Finite B = 4. 

^j^- = 100. In practice the value will be either much larger or much smaller than this 
number. In the plot, we still show the "misfitted" poles for reference. Overall, we see that 
the agreement of is very good at finite /, and the estimates from the spectral curves 
are much more reliable in the sense that there are no poles that are "not found" or have 
large displacements. In fact, even if there is some disagreement, the spacing Dq between 
the poles is more accurate than the overall shift of the poles. While there may be larger 
disagreements for the first pole and at small values of /, we should note that in those 
cases the direct fitting also fails frequently. There is no clear trend for the dependence 
on p and B. Typically the approximation is slightly worse for very large values and very 
small values, because in the former case the decay is more rapid, causing the resonances 
to be more asymmetric. In the latter case, we also see only the first few resonances and 
are not in the slowly decaying oscillatory regime - this time because of the larger spacing 
between resonances. We will not demonstrate this in the plots in fig. 30, simply because 
7o/^o is approximately constant when varying B or p at fixed /, and hence this is difficult 
to display in a clear fashion. The agreement in the imaginary direction is slightly worse, 
which is not unexpected because it depends more on the conjectured form of the residues 
for the Ansatz. 

Next, let us look at the spectrum of the resonances as a function of the resonance level 
n. In figure 31, we show the real part of the poles u n for various choices of p and B at 
M q = 0. Because it is important to distinguish between behaviors that are, for example, 
of the kind v/n(l + nj from strictly linear behaviors, we used the direct search for poles in 
order to obtain the highest accuracy. It turns out that the poles follow extremely closely 
a linear behavior v n = uDq + bv (or v n = (n — 1/2)Dq + bv in the case of a magnetic 
background field), with a small negative value for bv. The latter is easily explained from 
the behavior of the first pole that we found in the hydrodynamic regime. The constant 
separation implies that the spectrum of Landau levels is indeed the classical one for the 
finite mass case. This is slightly puzzling though, because one would have thought that if 
the collective excitations that carry the current have an effective mass, this mass would be 
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Figure 31: The real part v n of the location of the poles in the complex frequency plane as a 
function of the level number. The lines are linear fits. Left: Various values for B at / = 2. Right: 
Varying p at / = 4. 



finite. Then we should see a transition from the massive to the massless behavior at some 
frequency. Hence, either that mass must be large or frequency dependent, or some unusual 
mechanism gives rise to the Landau levels. 
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Figure 32: The real part ;>„ of the location of the poles in the complex frequency plane as a 
function of the level number at finite quark mass. Left: Varying m at B = 4 and / = 4. Right: 
Different values of p at m = 8 and f = 2. 



Looking at the case of finite quark mass in fig. 32, we find again no sign of a non- 
constant spacing between the poles, surprisingly even around the level of the quark mass. 
There is however a transition in the value of i>o around that region, which we will follow 
up on later. 

Finally looking at the behavior of the imaginary part of the poles, which reflects the 
inverse lifetime, in fig. 33, we find a small but significant deviation from the linear relation 
7n = n 7o + <$7 (and accordingly for the Landau levels). This small drift towards larger 7„ 
for small n can again be explained from the non-trivial behavior of the first (hydrodynamic) 
poles and from the fact that we are in a finite temperature background, which renders the 
first few resonances that are close to the temperature scale less stable. 
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Figure 33: The imaginary parameter in the location of the poles, uj n 
values of p at / = 4. Right: Varying mass m at / = 4 and B = 4. 



- *7«- Left: Varying 



Finally, we can look at the resonances. 
Before looking at Dq, let us study the shift 
5v in fig. 34. We find that for large /, there 
is a universal behavior bv oc B^ 1 , pT 1 with 
a proportionality constant that seems inde- 
pendent of /. A significantly different be- 
havior exists only for small values of / and 
large values of B or p. This may be sim- 
ply due to the worse fitting because in those 
cases we found only the first two poles, and 
the second one has already a very low am- 
plitude such that it is at the limit of what 
can be recognized as a resonance above the 
background. This overall behavior may be simply due to the fact that the first quasiparticle 
pole originates from the relaxation poles on the imaginary axis as we found in the previous 
section. 
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Figure 34: The overall shift in the poles 5v as a 
function of B or p for various values of /. 
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Figure 35: The separation Dq in the spectrum of "plasmons" and "Landau levels" as a function 
of B or p for various values of /. The plot on a right only shows the highest / curves. 
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Next, let us look at the value of v$ in fig. 35. There are two features to notice: Firstly 
the finite value of z>n at vanishing magnetic field or density, and secondly the non-trivial 
dependence on B or v - both of which are quite different from what we would have naively 
expected for plasmons and Landau levels. It turns out that there are two equivalent ways 
to interpret this situation. 
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Figure 36: vq presented as v 2 as a function of the square of the magnetization or chemical potential 
M 2 or jl 2 for various values of /. The plot on a right only shows the highest / curves. 

An interesting physical picture can be obtained by representing the data as z>q as a 
function of the square of the magnetization M or chemical potential fi that are computed 
in [31]. As we see in fig. 36 the result are perfectly straight lines, such that the resonances 



are given by . 



M 1 



or uoz 



Obviously, we could interpret the function 



a(f) as some kind of a mass scale that depends on / only, i.e. not on the width of the 
defect, but only its "topological" property. If we were - inappropriately - to look at the 
corresponding classical Schrodinger equation, we would see that then the magnetic and 
density perturbations are not independent, but mixes with some other potential. This is in 
contrast to the fact that the overall "amplitude" of the resonances, and hence the residue 
of the poles is at least for small fields proportional to the magnetic field. This is just what 
happens in the case of the classical Hall effect as discussed in section B.l. In fact, we can 
check a few values for uq, and compare them to the resonances at finite / due to a finite 
"width" of the defect, that were found in [28] to find that they are identical to the / — > 
limit of those reseonances. 

Inspired by this, we can try to check the relation between the resonances and the length 
scale given by the effective temperature that we observed in fig. 36. To do so, we can plot 

2. 



the ratio ^ — in fig. 37, where we see that it approaches within errors 



ML 

T eff 



This is not a big surprise, as there is an underlying exact relation between '^r- 



and 



fi (or M), that we can easily verify numerically and may in principle be able to derive 
analytically. 

Before concluding the study of the parameters of the resonances, let us look at the ratio 
3^ in fig. 38. Firstly, we notice the steps in the curves, that happen to coincide with data 
sets at which the highest pole in the sequence of resonances drops out of the fit because of 
its decreasing amplitude. This indicates the limitations in fitting the parameters (that we 
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Figure 38: The "inverse lifetime to mass ratio" 3^ for the magnetic and density resonances as a 
function of B or p for various values of /. 

expect to converge only asymptotically) accurately. Obviously, we could try to account for 
the non-linearity in j n as it was done in [28], but this has the downside that fitting with 
more parameters makes the result less reliable and may simply hide the limitations of the 
numerical result. Secondly, we see that within those limitations and even though ^ itself 
seems to approach constants in v or B, it seems to be best described as being proportional 
to Within the errors, it seems that the appropriate ratio is independent of / and 

becomes also approximately constant in v and B. 
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Figure 39: D for various values of the quark mass parameter rh at / = 4. Left: As a function of 
the magnetic field. Right: Vq as a function of the magnetization M 2 . 

Next, we can look at the location of the poles at a finite quark mass in fig. 39. If we 
look only at uq as a function of B, we see a surprise, as there is only an overall shift in 
the curves depending on the quark mass. On the other hand, if we plot v 2 as a function 
of the square of the magnetization M 2 , we see clearly the scale M ~ rh that separates 
the massive and massless regime. Below M ~ rh, i>o is suppressed with increasing mass 
and above M ~ rh, the behavior is similar to the massless case. If we compare z>o to 
the effective temperature, we find that is approximately constant around 2.04 . . . 2.06 
with no apparent systematic trend, so even in this case the effective temperature sets the 
appropriate scale for the quasiparticle energy spectrum. In principle, it would be interesting 



to study also smaller values of /, however at smaller / and large masses, it is not possible 
to find reliably at least the first two or three poles. 




Figure 40: Left: The difference in v$ between the case of the magnetic field and density at / = 4 
for various values of the quark mass as a function of p = B. Right: The scale T e ft/T at / = f and 
various values of the quark mass as a function of the net baryon density and magnetic field. 
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Figure 37: The ratio ^rjj as a function of B or 
p for various values of /. 



In principle, we expect that the results 
in the presence of the finite density and 
the magnetic field are different as we turn 
on the finite mass. Comparing v for these 
cases in the left fig. 40, we see that the 
difference is very small, even as m > B. 
We are uncertain as to whether these devi- 
ations are significant. From the behavior at 
fh = 8, it seems that there may be a small 
effect, which is suppressed in the quasipar- 
ticle regime. This question is resolved on 
the right in 40, where we plot T e ff/T at 
finite mass. We see that there is indeed a 



small difference between the dependence on p and B at finite quark mass. Upon close in- 
spection we also notice the separation between the regimes above and below the mass scale 
with the scalings oc B 1 / 2 or oc B, respectively, at least in the case of the magnetic 
field. It seems that there, the scale is B ~ fh, whereas for the density it is p ~ fh 2 . 

Because of the very good agreement between the quasiparticle spectrum and T e ff, 
let us have one more look at the mass-dependence of ^kf^ in fig. 41. Generically, the 
dependence on p and B is very similar at large values of / and differs very significantly 
at small values, hence our choice of plots. At small values of /, we see a very interesting 
behavior. It seems that at vanishing /, T e ff is essentially given by the mass, with only a 
subleading dependence on p at large masses, which actually reduces T e ff with increasing 
p. As we turn on /, this behavior turns over into a more common behavior, starting 
first at small / at large masses and small densities. This may be some transition from a 
purely 2-dimensional system to a system that extends also in the third dimension. What 
makes this behavior so surprising is that normally, both in the mass dependence and also 
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Figure 41: Left: -^r- as a function of the mass and the net baryon density for various (small) 
values of /. We choose the lower bound of the range 2 < p < 16 to avoid the phase transition at 
the critical mass and numerical problems near it. Right: as a function of the mass and the 

magnetic field for various values of /. 



at M q = 0, any dependence on / is subleading at large p. We should be careful with 
conclusions however, because we have not tested the dependence of the resonances on T e ff 
in this regime. Looking at the B dependence at large /, we find that the behavior is more 
generic, with a ~ S- scaling of T e jf atm>B and dependence approximately proportional 
to at m <C B. We also note that The transition between the two regimes becomes 
clearer at increasing /, i.e. as we widen the defect. Large / also suppress the B dependence 
ofT e// . 
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Figure 42: The "inverse lifetime to mass ratio" 
for / = 4 and various values of fa. 
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To conclude, let us also in this case look at the ratio ^ which we show in fig. 42. 
Again, we see that j£ is approximately constant in B, however it seems to depend on the 
quark mass. If we divide by the effective temperature, the dependence on the quark mass 
is removed, however there seems to be some dependence on B. 

Finally, we can look at how the poles split if we turn on the Hall effect. In fig. 43, we 
show the location of the poles for various values of /, p and B. In black, we show the poles 
in the absence of the Hall effect, and then we show the sequence of poles as we gradually 
turn on the "other" parameter. The most surprising result is that the total spacing Dq of 
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Figure 43: The splitting of the poles due to the Hall effect at various values of p, B and /. The 
black symbols indicate the sequence of poles in the absence of the Hall effect for some choice of / 
and B or p. With decreasing gray shade, we turn on the p or B, respectively, causing the original 
pole to split in two poles. 



the pairs of poles remains approximately unchanged, with quadratic dependence on the 
"smaller" background quantity 8vq(B) := \jv^ — vq(B = 0) 2 oc B that is within errors 
consistent with the behavior of the effective temperature. The splitting of the poles, i>A 
depends approximately linearly on the magnetic field or density that we turn on, however 
there seems to be no simple dependence of the proportionality coefficient on the obvious 
candidates such as the specific magnetic moment the density of states or the magnetic 
susceptibility. We show this in fig. 44 

The origin of the resonances on the gravity side is again straightforwardly explained 
in terms of quasinormal modes on the brane. Taking the same Ansatz A y = Aoe^ ds ^ as for 
the effective temperature, but now in terms of the variable s that we used in section 4.1.1 
and computing the equation of motion for £ as in (4.36), we obtain 

{y^GG tu G xy } J ^ - GG U G uu G xx G yy uj 2 = . (5.5) 

Since we just want to have a brief picture, we will only work to leading order, i.e. we use 
the approximate solution Co := iwv 7 GG a G uu G xx G yy and the pertubation ( = (o + e. This 
gives us the linearized equation of motion for e 



= e - 2eCo + Co =: e - ea(a) - f3(s) , 



(5.6) 




Figure 44: Left: The splitting of the energy of the quasiparticle poles due to the Hall ef- 
fect, i>a- Right: The shift in i> from turning on the second parameter, written as 8vq(B) := 

\Ji>l~ v {B = 0) 2 oc B (and p <-» B). 



again with the general solution e = e^o dsa ( s ) (eo + Jq ds eT fo dsa ( s "> [3(s)j . Essentially 
what happens now is that the resonances arise from the inhomogeneous term. If we imagine 
that the source term f3 were a delta function at some position so with amplitude £, then 
the contribution from this term would be e = ° dsa ( s ). Setting s = and taking the 
real part gives us the resonances we want as 5a ~ £ cos 2co V GG tt G uu G xx Gyy . 

To see where the resonances originate from in the geometry, we consider the high 
frequency limit Co 1 and remind ourselves that the integral of a periodic function vanishes, 
however the contribution to an integral of the type j F{x)e %wx from some region around xo 
will be of the order — . In our case, the frequency in the exponent is also not constant, 
but we can take care of that by a coordinate change s — ► j £o- After taking the derivative 
and changing back to s, we find that the term that gives us the contribution in the integral is 
ds d s . This expression already shows us straightforwardly that these resonances appear 
only in the presence of the background fields, since otherwise V ' —GG tt G uu G xx Gyy = 1. 
In fig. 45, we show this term for various choices of the magnetic field, / and the density. 
We see how this length scale arises, and we see also the structure that gives rise to the 
line splitting in the Hall effect. It also demonstrates how the different phase between the 
magnetic and plasma resonances arises, essentially through flipping the sign in this term. 
The shape of this contribution gives rise to the amplitude of the oscillations, and the 
positive and negative sections tell us that there is a higher suppression by Co. Note that at 

small u <C \J B, p we have s ~ u and ( ~ y/l + f 2 . Hence at small or 0(1) values for B 
and p, the spacing between resonances will be mostly controlled by the value of (, and not 
through s. The curves in the plot however control through their shape the nature and the 
stability of the resonances. 

Essentially, we can interpret this as scattering off a potential step, as it was dicussed 
in [28] and in a similar context als in [42]. 
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6. Discussion and Conclusions 
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In this paper, we used holographic tech- 
niques to investigate the transport proper- 
ties of certain defect CFT's. In particular, 
we studied matter on a (2 + l)-dimensional 
defect emersed in a heat bath of 3+1 di- 



mensional M = 4 SU(N C ) SYM theory. 
Compared to the analysis in [28], we added 
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hand it allowed us to compare the rich structure of our results to known phenomena in 
condensed matter physics, that are based either on generic considerations or on physics in 
the weakly coupled regime. We tried to distinguish between a) generic properties that seem 
to be independent of physical details of particular models, b) intuition that carries over 
from the weak coupling regime and c) properties that are specific to the strong coupling 
regime and allow us to build some new intuition that is generic in the strong coupling 
regime. 

Most of our analysis covers two distinct defect CFTs. The first was realized by em- 
bedding N { probe D5-branes in an AdSs x S 5 background and in this case, the system (at 
T = 0) preserves eight supersymmetries. The second system involves embedding vV f probe 
D7-branes in the AdSs x S 5 and the resulting defect CFT preserves no supersymmetries. 
As it contains only fermions at the massless level, it is of particular interest to condensed 
matter physics. In section 2.1, we looked at the gravitational setup that corresponds to 
turning on the various parameters in the field theory side. In both cases, the theory could 
be deformed by introducing a flux in the AdS part of the world volume corresponding to 
a finite net baryon number density and chemical potential in the dual CFT, a flux on the 
sphere corresponding to a shift N c — > N c + 5N C in the level of the gauge group on one 
side of the defect and a deformation of the sphere corresponding to a finite quark mass 
in the D5 case. In the D7 case, the internal flux was needed to stabilize the setup, and 
the operator dual to the deformation of the sphere seems to have non-integer conformal 
dimension as described in section 2.4. It turns out however that it is not quite sure in how 
far the D7 setup persists in the light of gravitational backreaction. 

Perhaps surprisingly, the transport properties of both defect CFT's were essentially 
identical in the massless case. Certainly, higher order effects might not be identical any- 
more. Furthermore it is curious that the only change is a Chern-Simons term with non- 
constant coupling at finite quark mass. However, we did not pursue this avenue at the 
massive level, since the results would be unreliable. 

We then went on to provide a few general results in section 3. In 3.1, we studied 
how the field theory outcome of the EM duality in gravity side changes in the case of the 
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Figure 45: Coefficient giving an estimate of the 
contribution to the resonances as described in the 



ground magnetic field, finite quark mass 
and net baryon density. This allowed us 
on the one hand to enlarge the class of the- 
ories that we are studying, and on the other 
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extra parameters. We found the very interesting result, that the transport properties are 
now related under a simultaneous exchange of the dimensionless magnetic field and den- 
sity and the transverse and longitudinal coordinate, while interchanging the 2-dimensional 
conductivity tensor with its inverse. From a condensed matter point of view, this related 
completely distinct parameters of the theory for a large class of 2+1 dimensional theories 
whose gravity dual obeys EM duality. Throughout section 5, we saw how this duality gets 
broken as we consider a finite mass. However, still, this breaking appears gradually, with 
the parameter M q /T. This duality seems to be part of the SL(2, Z) duality discussed in the 
field theory [44, 43], extended to finite frequencies and complex conductivities. The second 
generator may then be obtained from the theta term and "topological" Hall conductivity 
discussed in [28]. Certainly, it would be an interesting point to address in how far strongly 
coupled systems observerved in nature display such a duality and possess a self-dual point 
in phase space - in our system at the conformal point at vanishing net density and magnetic 
field. 

In section 4, we then discussed analytic results in several regimes. In the DC limit, 
we found the Drude conductivity at finite density, a magnetoresistance effect and a Hall 
conductivity. They can be parametrized under the Drude model, obviously giving a new 
description to the "Drude-parameters" in terms of the parameters of the theory, as the 
underlying microscopic physics is different. However, the overall scalings in the limit of 
large density and magnetic field are the familiar ones. Coincidentally, if we assume that 
the form of this DC conductivity has a high degree of generality, this could address the 
minimum quantum conductivity in Graphene [48] . There, it is known that at the neutrality 
point, which corresponds in our case to vanishing net baryon density, there is a minimum 
in the conductivity of e 2 /h per carrier type, and this has apparently been of significant 
interest in the community [47]. Under an appropriate translation of the parameters, this 
is precisely what we observe in our case, as we also find, an increasing conductivity as 
we move away from the neutrality point, which goes beyond the results in [10]. However, 
for example from the observed magnetic resonances, our defects seem to be quite differ- 
ent from the chiral nature of graphene [41]. We then studied the small- frequency limit, 
in which we reproduce the existence of a Drude peak and the minimum at lo = in the 
case of the magnetoresitance. We could also identify a relaxation time, and accurately 
reproduce a relation between the frequency dependence of the diagonal conductivity and 
the hall conductivity, in the limit of large densities. In general, the structure of the fre- 
quency dependence resembled the generic prediction from the Drude model. We also found 
however, that the behavior of the specific values of parameters that depend on a particular 
model changes compared to the Drude model - for example the relaxation time receives a 
dependence on the magnetic field. 

When we compared the relaxation time to the location of the purely dissipative poles 
in the hydrodynamic regime that we studied numerically in section 5.2, we found that there 
was a disagreement with the dominant "relaxation pole" only at small densities. This can be 
explained from the remarkable constant DC conductivity from the EM duality at vanishing 
density found in [28, 10], that is obviously not considered in the Drude model. In general we 
observed, in this regime and elsewhere, that there is always a total finite quark density in 
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thermal equilibrium (which we can't control) that influences the transport properties even 
at vanishing net density. Beyond the dominant relaxation pole, that showed an unusual 
dependence on the density, which we attributed to the strong coupling, there was a second 
relaxation pole, that shows a more classical behavior, but has no significant contribution 
to the charge transport. We were able to reproduce the transition to the quasiparticle 
regime that was found in [28], and also found another transition in which the relaxation 
poles merge and turn into the first Landau level at some critical magnetic field, rather 
than the diffusion pole merging with the first relaxation pole as in [28]. This transition 
can be attributed to strong coupling and is absent in free particle models. Overall, there 
is an interesting interplay between those relaxation poles and the diffusion pole as we tune 
the parameters, and there is a common theme that before poles move to large (imaginary) 
frequencies or leave their regime of validity, they either merge into quasiparticle poles or 
have decaying residue. 

The diffusion pole also has an unusual dependence on the density, which can be mo- 
tivated from the strong coupling properties. We verified that the numerically obtained 
diffusion constant agrees with the one obtained analytically from the membrane paradigm. 
We also computed the permittivity. This gives us what we called the "relative" permit- 
tivity, that depends on the mass, magnetic field, internal flux and density. The diffusion 
constant and permittivity reproduce from the diffusion behavior, i.e. from the Einstein 
relation, the precise value of DC conductivity. 

In the opposite regime, i.e. in the low temperature limit, our analytical approximations 
were concentrated at the exponentially suppressed conductivity at small frequencies and at 
large wavenumbers, q~^> T. There, we again found the conduction threshold at to = q and 
extended the result of the exponentially suppressed diagonal conductivity that was found 
in [28]. Now, the "effective temperature" that controls the "Boltzmann factor" however 
also depends on the other parameters of the theory. In particular, the density and magnetic 
field now raise that factor, i.e. reduce the exponential suppression. This is however not to 
be misunderstood as doping a semiconductor. At very large values of those parameters, 
of the order B, p S> q 2 , we were able to demonstrate that the conductivity turns into the 
DC result. Computing the Hall conductivity gave an interesting result as we obtained a 
finite value of the Hall conductivity even in the regime where the diagonal conductivity 
vanishes. This is however a common theme in condensed matter physics, for example in 
semiconductors at low temperatures or on Hall plateaus in the quantum Hall effect. 

The rest of the work in the quasiparticle regime was mostly numerics-based. In section 
5.1, we gave an overview over the spectral curves, where we found the appearance of what 
one could describe as the strong-coupling equivalent of Landau levels and plasmons, and 
line splitting in the Hall effect. We also noticed that increasing the mass or / makes the 
results approach the low temperature limit, as we expect from our results on the effective 
temperature. Furthermore, we looked at how the resonances from the isotropic regime carry 
over to the finite- wavenumber regime and connect to the resonances in the conformal case, 
which were discovered in [28]. In particular, we found that for the transverse correlator, 
the density resonances connect smoothly whereas the magnetic resonances connected less 
smoothly - which can be explained in terms of the localizing property of the magnetic 
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field. The other approach to the quasiparticle regime was to extract the location of the 
poles in the correlator in order to obtain the quasiparticle spectrum. We found that the 
poles are exactly equally spaced, indicating that the mechanism underlying the magnetic 
resonances is just a quantum harmonic oscillator as in the classical generation of Landau 
levels. The length scales corresponding to the spectrum can be explained in two manners: 
On the one hand, they are just given by approximately 2 times the inverse of the effective 
temperature, over essentially all the parameter range including the mass. On the other 
hand they can be related to the magnetization or chemical potential, which splits in the 
massive case into regimes below and above the quark mass. In this scenario, however, it 
seems that the Landau levels or plasmons are strongly coupled to the resonances over the 
"width" of the defect, as apparent from the minimum spacing of the resonances. This fits 
in nicely with the line splitting in the Hall effect, where each pole splits in two, indicating 
that there are overall two types of resonances in the system. Because of this and because 
of the unusual spacing, it seems that the magnetic and density resonances are not Landau 
levels or plasmons in the classical sense. In a qualitative description, we also discussed how 
the quasiparticles arise from quasinormal modes in the scattering off a step in a potential 
in the gravity side. 

Comparing our results to those obtained from field theory methods [18], we found 
that there were a few similarities as a resonance or threshold at u> = q is also generically 
obtained using field theory methods. Furthermore, our results can be expressed in terms 
of a universal function that depends on uj/T. However, it turns out that expressed in this 
way, this universal function depends on quantities like po/T 2 , and hence also depends on 
the temperature. Overall, it seems that using AdS/CFT, we could more straightforwardly 
obtain a very rich and complex behavior of this universal function. Also, it seems that the 
methods in [18] do not find quasiparticle resonances, that seem to be an integral part of 
the defect that we studied. 

For directions of future research, it would certainly be interesting gain a better inter- 
pretation of our results in terms of the microscopic theory beyond what we have attempted 
in this thesis. This reveals a big weakness in using the AdS/CFT correspondence, as it is 
in practice to some degree like performing an experiment; and the "microscopic" theory in 
AdS/CFT is not the field theory, but the gravitational configuration. Another interesting 
direction would be to study the problems of the D3-D7 defect more in detail, as having a 
purely fermionic system is very appealing, even if it is only in the sector of the fundamental 
representation on the defect and the 3+1 bulk is still SYM. It would be interesting to see 
what kind of effects may then appear in the "massive" case. 
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A. Components of the metric 



In this appendix, we write out the components of the metric in the most convenient com- 
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bination, including the implicit effective coupling in \J—G 
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At M q = 0, i.e. = 0, these simplify to: 
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B. Weak-Coupling Condensed Matter Physics 

Even though we are interested in the strong coupling regime which one expects to be 
quite different from the free electron gas picture, some intuition and generic properties 
can be learned from this very straightforward limit and it can serve as a phenomenological 
description. Also, weak coupling is the only reference regime over which we have good 
control and where there are readily available textbook-type results. Hence, we remind the 
reader of the very basic model, which can be found in standard textbooks [39] . 



B.l Metals 

In the Drude model, we assume a gas of non-interacting charge carriers with finite (effec- 
tive) mass m e ff and charge e, which we will write out explicitly. Eventually, it will turn 
out, however, that using the coefficients that do not involve m e ff is a suitable parametriza- 
tion also in the relativistic case. To obtain the conductivity, one then considers a small 
electromagnetic background field, to which the charge carriers are coupled classically via 
the Lorentz force m e ff(dtp— T^p) = zE + ev x B. As the charge carriers are massive, 
they have a finite net velocity v, which is assumed to be neutralized on the time scale of a 
relaxation time r. Classically, one has then a mean velocity v = m e ^ E =: fj,E, where we 
defined the charge carrier mobility [i. The charge carrier mobility is related to the Diffusion 
constant by the Einstein relation 

D = — . (B.l) 
e 

Further, for our massive case, the magnetic field can be rewritten in terms of the cyclotron 
frequency as uj c = = \*B /r. Now, let us consider two species of charge carriers with 
equal mass and relaxation time, but opposite charge ±e, such that we have a total density 
of charge carriers n = n + + n_ and a net charge density An = n + — n_ . This is relevant in 
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our case, since even at vanishing net baryon number density, po = 0, at finite temperature 
T > M q we will always have a finite total baryon density. 

To obtain the conductivity, one then assumes an oscillatory electric field E = Eoe~ lujt 
and current j = joe~ lwt , but constant magnetic field B and obtains the diagonal conduc- 
tivity 



II en^(l - iurr) en/x / , 1 - uj 2 t 2 2 2 1 - 3u; 2 r 2 

" (1-^t) 2 +^ 2 t 2 l + o; 2 r 2 V WT l + ^r 2 WT (1 + w 2 t 2 ) 2 

(B.2) 

Taking into account the positive and negative charges, the Hall conductivity becomes 
eAn^ c r _ eAn^r f 1 + _ 2zog_ _ ^ _3-_c^ + ^3 

(B.3) 



(1 - ioJT) 2 + w 2 r 2 1 + co 2 r 2 V 1 + w c 2 r 2 (1 + oj 2 t 2 ) 2 

The dissipative part of the conductivity is then 



II 1 +0J 2 T 2 +U) 2 T 2 

Re a = C¥ (TT^WW and (B - 4) 

o -L A 1 - W 2 T 2 + W 2 T 2 

Re a = eAn/xw c r— — — . B.5 

(1 + w^r^ — cj^t^)" 1 + 4lo^t z 

The DC conductivity at B = is commonly referred to as the Drude conductivity, and 
at small frequencies it is also called the Drude peak, due to the small value of r in metals 
at room temperature. Similarly, the fact that the DC conductivity is suppressed at finite 
magnetic fields is referred to as the magnetoresistance effect. 

It is also interesting to notice that at the magnetic resonance around cj c , the Hall 
conductivity changes sign, and this turns into a pole at large r, i.e. in practice at very 
small temperatures or in very "clean" semiconductors. 

We can also observe a few generic properties of the frequency dependence. For example, 
at B = 0, = ^ and, provided the relaxation time is independent of the magnetic 

field, the behavior at large frequencies is = — An interesting relation is also 

"f^fx = | a t vanishing magnetic field and 3 at large magnetic fields. Also, the magnetic 

field at which changes sign is t 2 oj 2 = \ for cr" and 3 for . The fact that the quadratic 
term changes sign implies that the Drude peak moves away from the real axis and becomes 
a magnetic resonance. 

We already see a limitation of the free electron gas picture, because - as we will see in 
section B.3 - the quantum mechanical treatment implies that the first magnetic resonance 
is at \oj c and taking into account a finite coupling implies that there will be plasma density 
resonances. 

B.2 Semiconductors 

Semiconductors are somewhat less generic than the Drude model of conductivity, i.e. we 
must assume that we are dealing with fermions, but the discussion obviously carries over to 
any system of charge carriers with an excitation gap. In a semiconductor, we assume that 
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the valence and conduction bands are separated, E v < E c , where the valence band is, as 
the name says, filled such that the chemical potential lies between the bands E v < \x < E c . 

At small temperatures, the dissipative conductivity will be dominated by n, such that 
we are most interested in obtaining the density of conduction "electrons" n c and valence 
"holes" n v assuming some density of states g c and g v : 

n c = £d Egc (E) e{E J /T + i (B.6) 



and 



^ = / EV dE 9v (E) (l - e(E J /T + 1 ) = J E JEg v (E) e{ ^ /T + i . (B.7) 



Assuming that we are dealing with low temperatures E c — fj, 3> T and fj, — E v 3> T, we can 
re-write this as 

n c (T) = N c (T)e~^-^ T , n v (T) = N v {T)e'^- E ^ )/T (B.8) 
where we defined the edge densities of states as 

[•oo pE v 

N C (T) := dEg c (E)e-( E - E ^' T , N V (T) :=/ dE g v (E)e~ (E ^/ T . (B.9) 

J E c J oo 



Using those definitions, we can combine the equations (B.8) to write down the "law of 
mass action" 

n c n, = N c N v e~ {E ^ E ^' T = n,(T) 2 (B.10) 

which determines the charge carrier density of an intrinsic (undoped) semiconductor, rij = 
n^, = n c . This small charge carrier density implies (together with the purity of the crystal) 
that r may diverge at small temperatures, as it is usually (in a metal or metallic phase 
with high charge carrier density) dominated by the charge carriers and their thermal motion 
(independent of the purity). Hence semiconductors may show, for example, a finite Hall 
conductivity even if the material is effectively an insulator. 

The chemical potential can then be obtained from the edge densities of state 

„_^ + . riB * (R11) 

One can see that in practice in an intrinsic semiconductor at E c — E v T, the edge 
density of states has only a very small influence on the chemical potential, and hence on 
the conduction threshold and the suppression of the conductivity. 

B.3 Resonances 

In the "optical" regime at larger frequencies, we are interested in quasiparticle resonances, 
in particular in Landau levels and plasmons. In the case of massive charged particles, 
the derivation of Landau levels is straightforward. Assuming coordinates {t, x, y, z} with 
the magnetic field in the z direction, the gauge potential can be written as A = B x dy. 
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Substituting this into the Schrodinger equation gives rise to a quantum harmonic oscillator 
with a frequency uj c = with the solution 

E„ = E%l + u c (n-±), neZ + . (B.12) 

(z) 

Obviously, our charge carries are constrained to the z = plane, so E y ki ^ = 0. Using the 
same naive strategy for a scalar in the Klein-Gordon equation, we find that 

E 2 n = p 2 z + m 2 + u 2 M (2n - 1) , uj 2 m = tB , (B.13) 

putting the Landau levels in the massless case at E n = ±cjmv / 2^ — 1- In either case, one 
can apply a simple argument by assuming a finite sample size to derive the density of states 
(per unit area) of 

N ± a ^ = tB . (B.14) 
A 

In a similar fashion, Landau levels can also be obtained from the Dirac equation (see [41] 
and references therein): For massless chiral fermions e.g in graphene, the result is [41] 



E n = ±lu m\v 'f|v \ n \ ) n£Z (B.15) 

where vf is the Fermi speed (vf = 1 for a "real" relativistic system) and for a system of 
chiral fermions with finite effective mass, as in multi-layered graphene, [41] 

E n = ±u cy /\n\(\n\ + 1) , neZ . (B.16) 

The zero-energy, field-independent Landau level is unique to chiral fermions, which makes 
them very easy to identify [41]. 

The other quasi-particles that we are interested in are plasmons. These are collective 
excitations resulting from density perturbations of the (electron) gas. They can be derived 
in several ways. Classically it can be derived simply from the continuity equation V-j = dtp 
and from Gauss' law V ■ E = 4np. Using a^E = j, we can arrive with an equation for uj, 
47rn / LiT~ 1 = lo(lu + ir _1 ). Another classical derivation is to assume a neutral gas of positive 
and negative charges in which the charges are displaced in, say the x direction, leaving 
two strips of density ±An with a width d that obeys the classical equation of motion 
d 2 d = — 47re J fc A "^ . Finally, a proper derivation is based on computing the Green's function 
in a gas of weakly interacting fermions (of equal charge in magnitude and sign) with the 
Coulomb potential. A very instructive derivation can be found in [40] and yields 

. p = ^ = ^ . (B.17) 

where /cp is the momentum corresponding to the Fermi energy. The Green's function is at 
small momenta \q\ <^ \u\ 

D = % — 2 . (B.18) 

1 (uJ + iT- 1 ) 2 
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The derivation follows with only small changes also in the massless case, where we find to p = 
kp8ir/3 . Plasmons are usually (i.e. in weakly coupled systems) observed through optical 
scattering, where one can observe spectra from multiple plasmon excitations. Studying 
plasmons and surface plasmons in various materials an using them for photonic devices 
seems to be a very active field of research. 
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